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ABSTRACT 


A  mathematical  model  was  formulated  to  describe  the 
injection  of  saturated  steam  down  oilwell  tubing  under  constant 
inlet  conditions. 

The  system  was  divided  into  three  parts,  consisting 
of  the  fluid,  the  wellbore  and  the  formation.  Each  part  was 
considered  to  be  a  separate  system,  with  the  heat  flux  across 
the  boundaries  acting  as  the  coupling  parameter.  The  wellbore 
consisted  of  a  tubing  string  enclosed  by  either  one  or  two 
casing  strings.  The  heat  transfer  in  the  wellbore  was  con¬ 
sidered  steady  state.  The  heat  transfer  in  the  formation  was 
treated  as  unsteady  state  radial  conduction  and  the  fluid  flow 
was  described  by  a  modified  two  phase  flow  correlation. 

The  analysis  resulted  in  three  equations,  two  of  them 
implicit,  which  were  subsequently  solved  simultaneously  on 
the  7040  computer.  The  solution  did  not  present  any  con¬ 
vergence  problems.  Although  a  complete  verification  of  the 
mathematical  model  was  not  possible,  the  calculated  tempera¬ 
ture  profiles  were  similar  in  shape  to  an  observed  temperature 
profile . 
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INTRODUCTION 


With  the  cost  of  finding  and  developing  high  grade 
crude  reserves  continually  increasing,  more  attention  is  being 
focused  on  the  development  of  high  viscosity  or  "heavy"  oil 
reserves.  Since  the  oil  viscosity  is  highly  dependent  on 
temperature,  a  small  increase  in  reservoir  temperature  de¬ 
creases  the  oil  viscosity  markedly,  facilitating  increased 
production  rates  and  ultimate  recoveries.  One  of  the  more 
popular  means  of  increasing  the  reservoir  temperature  is 
steam  injection.  To  evaluate  the  feasibility  of  a  thermal 
project,  a  reasonably  accurate  estimate  of  the  heat  losses 
from  the  wellbore,  sandface  pressure  and  sandface  quality  of 
the  steam  would  be  desirable,  if  not  necessary.  Since  the 
injection  times  are  generally  quite  small,  steady  state  solu¬ 
tions  would  not  be  very  representative  and  unsteady  state 
solutions  are  required. 

The  wellbore  heat  transfer  coefficient  contains 
radiation  and  natural  convection  terms,  both  of  which  are 
temperature  dependent.  The  nature  of  the  radiation  and 
natural  convection  contributions  present  major  difficulties 
in  deriving  an  unsteady  state  equation  describing  the  heat 
transfer  in  the  wellbore,  and  it  appears  that  a  numerical 
solution  is  the  best  that  could  be  hoped  for.  To  avoid 
lengthy,  time  consuming  numerical  techniques,  the  problem  was 
divided  into  parts;  whereby,  wellbore  heat  transfer  was  con- 
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sidered  steady  state  and  formation  heat  transfer  was  treated 
as  unsteady  state  radial  conduction.  The  fluid  had  to  be 
considered  since  flowing  friction  affects  the  fluid  tempera¬ 
ture,  which  in  turn  affects  the  heat  losses.  The  three  parts 
were  then  solved  simultaneously,  with  the  heat  flux  across 
the  various  boundaries  serving  as  a  common  boundary  condition. 
(Heat  lost  by  the  fluid  =  heat  transferred  through  the  well¬ 
bore  =  heat  entering  the  formation.) 

The  heat  losses  result  in  an  enthalpy  change  in  the 
fluid  causing  a  temperature  drop  in  a  single  phase  fluid  or 
a  quality  change  in  a  two  phase  steam.  Other  investigators 
(11,  13,  15)  have  treated  the  single  phase  case.  Satter(13) 
considered  the  two  phase  case  neglecting  friction  and  kinetic 
energy  effects,  thus  obtaining  a  constant  steam  temperature 
in  the  tubing.  Under  certain  conditions  this  is  a  poor 
assumption,  as  was  substantiated  with  field  data  (3)  (see 
Figure  4) ,  where  a  temperature  drop  of  44°F  was  recorded  for 
a  1600  ft.  we 1 1 . 

This  investigation  was  limited  to  an  attempt  to  treat 
the  injection  of  a  saturated  steam  at  constant  inlet  conditions 
(quality,  flowrate,  temperature)  as  fully  as  possible.  The 
resulting  generalized  computer  program  is  subject  to  the 
assumptions  and  approximations,  which  are  believed  to  be 
realistic,  made  in  obtaining  the  mathematical  model.  The 
accuracy  of  the  results  is  largely  dependent  on  the  accuracy 
of  the  empirical  correlations  employed  in  the  solution. 


;  r  j  a,'-  r  ,  ", 
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LITERATURE  REVIEW 


Two  authors  who  have  concerned  themselves  with  the 
specific  problem  of  heat  losses  from  the  wellbore  during  hot 
fluid  injection  are  Ramey (11)  and  Satter(13).  Since  their 
work  formed  the  basis  of  this  investigation,  it  is  treated 
thoroughly  in  this  review. 

Ramey  was  primarily  concerned  with  the  temperature 
profile  during  single  phase  fluid  injection.  Application  of 
his  model  to  saturated  steam  leads  to  a  simple  relationship 
subject  to  considerable  error.  Ramey's  basic  assumption  may 
be  listed  as  follows: 

1)  the  physical  properties  of  the  fluid  and  the  for¬ 
mation  are  independent  of  depth  and  temperature. 

2)  the  heat  transfer  in  the  wellbore  is  rapid  compared 
to  the  rate  of  heat  transfer  in  the  earth,  and  may  be  con¬ 
sidered  steady  state. 

3)  the  overall  heat  transfer  coefficient  is  independent 
of  depth,  and  may  be  determined  prior  to  a  temperature  cal¬ 
culation  . 

4)  the  frictional  losses  and  the  kinetic  energy 
effects  are  negligible. 

With  the  application  of  energy  balances,  simplified 
by  the  stated  assumptions,  Ramey  (1.1)  developed  equations 
describing  the  temperature  of  the  injected  fluid  as  a 
function  of  depth.  For  the  liquid  he  obtained: 


' 
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T  =  a  +  bz  -  bA  +  (T  +  bA  -  a)  e 
s  s 

o 


•z/A 


and  for  the  gas  obtained: 


(1) 
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a  +  bz  ~  A  (b  + 


-)  + 


JC 


T  -  a 
sQ 


where 


+  A  (b  + 


JC 
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■z/A 


=  a  +  bz 


(2) 


(3) 


A  = 


W  .  C  (k„  +  r  U  f  (t)  ) 
st  p  E  c  o 


2  it  r  U  k^ 
c  o  E 


(4) 


The  time  function  f(t)  is  a  function  introduced  to  make 
equation  9  valid.  Ramey  presented  a  graph  of  f(t)  versus 
dimensionless  time  (tD) ,  with  UT  as  a  parameter.  (Refer  to 
Appendix  I-C  for  the  derivation  of  f(t)) 

The  equations  were  developed  applying  an  energy 
balance  on  an  element  of  fluid  of  length  dz  at  a  depth  of  z 
ft.  in  the  tubing,  and  equating  the  heat  lost  by  the  fluid  to 
the  heat  transferred  through  the  wellbore  and  the  heat  enter¬ 
ing  the  formation. 

Thus 


AQ1  =  AQ2  =  AQ3 


(5) 


t  Refer  to  the  "Nomenclature"  for  the  definition  of  terms. 
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where 

AQ^  =  the  heat  lost  by  the  fluid  (Btu/hr) 

AQp  =  the  heat  transferred  through  the  wellbore  (Btu/hr) 
AQ3  =  the  heat  entering  the  formation  (Btu/hr) 

For  an  incompressible  liquid  flowing  through  a  constant 
diameter  tubing,  the  energy  balance  reduces  to: 

g  dz 

W  (dH  -  — - )  =  -  AQ-,  (6) 

s ,  ,  1 

t  q  J 

Jc 

Assuming  that  Cv  -  and  that  the  friction  losses  are  negligible, 
application  of  thermodynamics  yields: 


Also  have 


C 

P 


(7) 


AQ0  =  2-n  r  U  (T  -  T  )  dz  (8) 

2  c  o  s  c 

and 

2tt  k  (T  -  T  ) 

AQ  = - £ - 5 - (9) 

f(t) 


Equating  Aq^  and  AQ^,  solving  for  T  ,  substituting  the 

expression  for  T  into  equation  8  and  subsequently  equating 

c 

AQ^  and  AQ2  yields : 

3  T  T  T„ 

- -  +  —  -  Jl  =  0  (10) 

3  z  A  A 

where  A  is  defined  by  equation  4  and  /  0.  The  solution  of 
equation  10,  subject  to  the  boundary  condition  T  ( 0 , t )  =  T 

b  b 

o 


f 
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results  in  equation  1. 


For  an  ideal  gas,  dH 


C  dT  ,  thus : 
p  s 


g  dz 


g  J 
^  c 


(11) 


Substituting  equation  11  for  equation  7,  and  following  the 
same  procedure  as  outlined  for  the  liquid  case,  results  in  the 
derivation  of  equation  2. 

The  rate  of  heat  loss  per  foot  of  wellbore  may  be 
expressed  by: 


2"rc  kE  (TS  -  V 

r  U  f (t )  + 
co  E 


(12) 


dz 


If  T  and  U  are  known,  as  would  be  the  case  for  the  injection 
so  J 

of  a  saturated  steam  with  negligible  friction  losses,  the 
calculation  of  q^  at  any  depth  is  straight  forward. 


Satter(13)  improved  Ramey's (11)  solution  by  circum¬ 


venting  the  restrictions  of  a  constant  wellbore  heat  transfer 
coefficient  and  temperature  independent  fluid  properties.  By 
limiting  the  integration  of  equation  9  to  small  increments 
in  depth,  Satter  could  choose  increments  over  which  the  well¬ 
bore  heat  transfer  coefficient  and  fluid  properties  are 
relatively  constant.  The  overall  wellbore  heat  transfer  co¬ 
efficient  was  represented  by: 
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where 

hc  =  convection  coefficient  for  the  annulus 

=  F  (T  ,  T  ) 
a  s  c 

h^  =  radiation  coefficient  for  the  annulus 

=  F,  (T  ,  T  ) 
b  s  c 

For  single  phase  flow,  Uq  is  clearly  dependent  on  the 
fluid  temperature,  while  the  fluid  temperature  is  dependent  on 
(see  equation  1  and  2) .  A  trial  and  error  solution  is 
required  to  obtain  a  value  for  the  fluid  temperature  and  the 
wellbore  heat  transfer  coefficient  simultaneously.  Satter 
solved  for  the  temperature  profile  by  straight  iteration, 
whereby  the  first  trial  temperature  is  assumed,  then  a  tempera¬ 
ture  is  calculated.  If  the  assumed  and  calculated  temperatures 
do  not  agree  to  within  a  specified  error,  the  calculated  value 
becomes  the  new  trial  value. 

For  saturated  steam  Satter  obtained  relationships  for 
the  quality  versus  depth  and  percent  heat  loss  versus  depth. 

From  equation  6  obtain 


AQ 


1 


=  dH 


W 


sic 


g  dz 


g  J 
J  c 


(14) 


For  a  saturated  steam  at  a  constant  temperature 

dH  =  H  •  dy  (15) 

v  1 

The  substitution  of  equation  15  and  equation  12  into  equation 
14  yields  : 
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B'  + 


T  -  (a  +  bz) 
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0 


(16) 


A' 


where 


W  ,  H  (r  U  f(t)  +  k„) 
st  v  c  o  E 


A' 


(17) 


2  tt  r  U 

c  E  o 
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B  ' 


(18) 


g  H  J 
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T 
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a  +  bz 


The  integration  of  equation  16  from  z  ~  Az  to  z,  yields  the 
equation  for  the  steam  quality. 


bAz2 


y  (z , t )  =  y(z  -  Az,t)  + - 

2A' 

[a  '  B  '  +  a  +  b  ( z  -  A  Z )  -  T  ]AZ 
+  - - — — - — 

A' 

Setter (13)  defined  the  percent  heat  loss  as  the  heat 
loss  divided  by  the  usable  heat  content  of  the  inlet  steam. 
The  usable  heat  content  was  defined  as  the  enthalpy  of  the 
inlet  steam  minus  the  enthalpy  of  the  liquid  at  the  surface 
temperature  of  the  earth.  For  the  injection  of  a  saturated 
steam,  the  mathematical  equation  would  be: 


Percent  heat  loss 


(20) 


Q 


Squier  et  al(15)  derived  an  analytic  solution  for  the 
temperature  profile  of  a  water  injection  well.  As  they  assumed 
an  infinite  heat  transfer  coefficient  between  the  fluid  and 
the  formation,  their  solution  would  be  limited  to  problems 
involving  injection  down  a  casing. 
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THEORY 

The  equations  stated  in  this  section  are  derived  in 
Appendix  I,  Sections  A,  B  and  C,  and  are  subject  to  the  fol¬ 
lowing  basic  assumptions: 

1)  saturated  steam  is  injected  through  the  tubing  at 
constant  inlet  conditions  (quality,  flowrate  and  temperature) 

2)  heat  transfer  in  the  wellbore  is  steady  state, 
whereas  heat  transfer  in  the  formation  is  unsteady  state 
radial  conduction* 

3)  physical  and  thermal  properties  of  the  formation 
are  independent  of  depth  and  temperature. 

4)  the  annulus  contains  air  at  a  constant  pressure. 

The  basic  difference  between  this  model  and  those  pro¬ 
posed  by  Ramey (11)  and  Satter(13)  is  the  inclusion  of  friction 
and  kinetic  energy  effects.  To  obtain  the  mathematical  model, 
the  physical  system  was  divided  into  three  systems,  resulting 
in  three  separate  equations.  The  three  systems  were: 

a)  the  fluid 

b)  the  wellbore 

c)  the  formation 

The  heat  flux  across  the  system  boundaries  served  as  the  para¬ 
meter  interrelating  the  three  systems.  The  equations  were 
derived  assuming  that  the  systems  were  uncoupled,  whereby  all 
the  variables  pertaining  to  another  system  were  assumed  to  be 
known.  The  final  solution  was  the  result  of  successive 
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FIG.  1.  DIAGRAM  OF  THE  PHYSICAL  SYSTEM 


Steam  at  constant  quality,  rate  and  pressure 


Well  Bore 
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iterations  until  all  systems  were  satisfied  simultaneously, 
a)  The  Fluid 

The  energy  balance  for  the  steam  yields: 


dH 


dz 

—  + 
J 


udu 


Assuming  the  casing  temperature  to  be  a  linear  function 
of  depth  over  the  interval  Az,  an  equation  describing  the 
steam  quality  at  the  bottom  of  the  interval  may  be  derived 
(see  Appendix  I-A) . 
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Y' 
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(Tc  "  TF  )  +  (Tc  "  TE  } 

-  1  1  2 


2  -1 


(A4) 


W  .  0  1 

f  =  ( - £t -  )  2  - 

SGOOTrr-,^  2 g  J 
1  Jc 


The  force  balance  relates  the  friction  losses  of  the 
steam  to  a  pressure  drop.  The  equation  describing  the  pressure 
drop  for  two  phase,  vertical  steam  flow  down  a  circular  pipe 
was  taken  from  Perry(lO). 


It  may  be  written  as: 
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.1- 
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fGPGUG 


dz 


gcrl 


(A8) 


where 

=  two  phase  friction  factor  =  f(REG,Y)  (A9) 

(The  graphical  relationship,  as  presented  in  Perry  (10)  is 
included  in  Appendix  III-A  as  Figure  14) 
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Equation  A8  may  be  rewritten  to  yield 


dP 


v. 


X  f  y  dz 


where 


W 


s  t 


(12-  3600)  2  g  7r2ri  5 
^c  1 

To  eliminate  the  need  for  interpolation  and  extrapolation 
routines  necessary  if  the  graph  for  fn  were  to  be  applied 
directly,  simplifying  assumptions  were  made.  The  liquid 
density  and  vapor  viscosity  were  assumed  constant,  and  f_  was 
plotted  versus  quality  for  various  inlet  conditions  (See 
Figure  15  and  Figure  17  in  Appendix  III-A) .  A  straight  line 
fit  was  generally  sufficiently  accurate  over  the  anticipated 
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quality  ranges.  At  high  injection  rates  the  coefficients  for 
the  line  were  relatively  independent  of  superficial  velocity 
(see  Figure  17).  Thus: 

f„  =  a  +  b  y  (A.14) 

G  Y  y 

Keenan  and  Keyes  (8)  present  mathematical  relationships 
between  vapor  pressure  and  temperature,  and  vapor  volume  and 
temperature  for  saturated  steam. 

P  =  4>(T  )  (A15 ) 

s 


V_  =  ip(  T  ,  4>(T  )) 

vjr  S  S 


(A16) 


employing  equations  All,  A14,  A15  and  A16,  it  was  possible 
to  derive: 


where 


X  (a  +  b  y)  v 

y  y^  J 


r  (T  ) 

- — -  dT 

^  (T  ,  <KT  )) 
s  s 


( Al  8 ) 


(A19) 
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The  substitution  of  equation  A7  into  equation  A18  results  in 

one  implicit  equation  in  z?  and  T  „  Upon  specification  of 

z  s2 

T  ,  z„  may  be  obtained  by  trial  and  error,  using  Newton's 
s2  ^ 

method  for  convergence. 
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b)  The  Wellbore 


Assuming  steady  state  heat  transfer  in  the  wellbore, 
the  heat  transferred  through  any  part  of  the  wellbore  must  be 
equal.  Based  on  this  criteria,  an  overall  heat  transfer 


t 

coefficient  based  on  the  outside  casing  area  may  be  defined. 
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where  k  is  defined  by  equation  B5  in  Appendix  I-B.  Equation 
c 

B3  may  be  rewritten  as  : 
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t  Although  only  the  single  casing  case  is  treated  in  this 

section,  single  and  double  casing  relationships  are  derived 
in  Appendix  I-E 
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Equating  the  heat  transferred  through  the  wellbore  to  the  heat 
entering  the  formation  and  solving  for  the  outside  casing 
temperature  yields: 


T 


kF  (Ts  "  V 

T  -  T  +  — — - - - =  0 

kE  +  U  rcf(t) 


(B7) 


Assuming  ~  A^_  and  +  A^,  where  A^  and  A^  are 


t 


independent  of  T  ;  the  substitution  of  equation  B4  into  equation 

c 

B7  yields  one  implicit  equation  in  T  .  The  casing  temperature 

(T  )  may  be  determined  by  trial  and  error,  employing  Newton’s 
c 

method  for  convergence.,  Thus: 


(n+1) 


T 


(n) 


F 


(n) 

T 


dFm  ,  » 
/ _ T.  (n) 


■dT 


After  a  value  has  been  obtained  for  T  which  satisfies 

c 


equation  B7,  Uq  may  be  calculated  using  equation  B4. 


c)  The  Formation 


Assuming  a  time  independent  steam  temperature,  the 
partial  differential  equation  describing  radial  heat  con¬ 
duction  in  cylindrical  coordinates  may  be  written  as: 


where 


3  2T*  1  3T*  1  3T* 

3r2  r  3r  a  3t 


T  * 


T 

s 


(C4 ) 
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Carslaw  and  Jaeger (4)  provide  an  analytical  solution  for  the 
partial  differential  equation,  subject  to  the  following  con¬ 
ditions  : 

Initial  condition:  T*(r  ^  r  ,0)  =  T  -  T  =  T  * 

C  s  o 

Boundary  conditions:  1)  T*  is  finite  at  r  =  00 


2) 


3T* 


9  r 


r=r 


T  * 


r=r 

c 


c 

Assumption:  Uq  is  not  a  function  of  time. 

Since  only  the  heat  flux  at  the  wellbore  interface  is  re¬ 
quired,  the  solution  in  Carslaw  and  Jaeger  may  be  simplified 


to 

yield : 

9 

r  -V2 

e 
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From  the  definition  of  f  (t) ,  the  following  relationship  may  be 
derived : 


(C5 ) 


(C7 

dw 

0) 


f  (t) 


T 

s 


r=r 


(CIO) 


The  evaluation  of  the  integral  in  equation  C7,  necessary  in 
the  calculation  of  f(t),  is  described  in  Appendix  I-C 

Summarizing  the  equations  developed  in  this  section, 
it  is  possible  to  describe  the  physical  system  with  three 
interrelated  equations. 


. 


. 
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1)  z  = 

Fl(Ts'  V  f  (t) '  z) 

(equation 

A18) 

2)  T 

c 

F 2 (Ts ,  f (t) ,  Tc,  z) 

(equation 

B7 ) 

3)  f(t) 

=  F3(Tc,  Tg,  t,  z) 

(equation 

CIO) 

As  there  are  five  variables  (T  ,  T  ,  f (t) ,  z  and  t)  and 

s  c 

only  three  equations,,  two  variables  may  be  specified.  Due  to 

the  nature  of  the  equations,  the  specification  of  T  and  t 

s 

would  lead  to  the  simplest  calculations.  After  all  variables 
have  been  determined,  other  parameters  such  as  quality  and 
heat  losses  may  readily  be  calculated. 


. 


2  9 


DISCUSSION 

a)  Discussion  of  the  Model 

The  basic  assumptions  leading  to  the  mathematical 
model  are  restrictive,  but  often  necessary.  Although  it  may 
be  possible  to  remove  them  by  a  completely  numerical  solution, 
the  complexity  of  the  problem  increases  by  such  an  extent 
that  any  gain  from  a  better  model  is  soon  lost  to  round-off 
errors.  It  is  believed  that  a  comprehensive  solution  circum¬ 
venting  the  assumption  of  steady  state  heat  transfer  is  dif¬ 
ficult,  if  not  impossible.  The  complicated  boundary  conditions, 
added  to  temperature  dependent  radiation  and  convection  terms, 
tend  to  rule  out  an  unsteady  state  analytic  solution.  The 
number  of  grid  points  for  a  numerical  solution  are  limited, 
as  radiation  and  natural  convection  make  it  difficult  to  sub¬ 
divide  the  annulus.  Comparing  the  heat  absorbed  by  the  well¬ 
bore  to  the  heat  transferred  through  the  wellbore,  it  is  con¬ 
ceivable  that  the  error  introduced  by  this  assumption  is 
negligible,  except  for  very  short  times. 

Constant  inlet  conditions  were  necessary  insofar  as 
they  guaranteed  a  relatively  constant  temperature  profile  for 
the  steam,  a  stipulation  imposed  by  the  analytic  solution  for 
f(t).  No  attempt  has  been  made  to  obtain  or  develop  an  analy¬ 
tic  solution  for  f(t)  incorporating  a  time  dependent  steam 
temperature,  but  it  is  unlikely  that  such  a  solution  exists. 


' 
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The  problem  was  limited  to  a  solution  for  saturated 
steam  injection,  as  existing  solutions ( 11 , 13 , 15 )  for  single 
phase  fluid  injection  are  believed  to  be  adequate.  For  single 
phase  fluids,  the  friction  loss  would  affect  the  heat  loss  to 
the  extent  that  pressure  would  affect  the  fluid  enthalpy. 

Since  the  pressure  drop  for  a  single  phase  fluid  is  consider¬ 
ably  lower  than  for  a  two-phase  fluid,  the  overall  effect  of 
friction  is  believed  to  be  minor. 

The  assumptions  of  constant  annular  pressure  and  depth 
independent  thermal  and  physical  properties  of  the  formation 
could  readily  be  removed  by  stratification  of  the  wellbore  and 
the  formation. 

The  additional  assumption  of  a  constant  steam  tempera¬ 
ture,  as  made  by  Ramey (11)  and  Satter(13)  eliminates  one 

equation  (A18)  and  one  variable  (T  ).  The  assumption,  result- 

s 

ing  from  the  neglection  of  friction  losses,  is  necessarily 
poor  for  flow  involving  two  phases  and  high  velocities.  Econo¬ 
mic  considerations  dictate  high  injection  rates,  resulting  in 
high  velocities , which  subsequently  cause  sizeable  friction 
losses.  Field  measurements ( 3)  verify  this  hypothesis  by  pre¬ 
senting  a  pressure  drop  of  100  psi  in  1600  ft.  of  tubing  (see 
Figure  4).  Neglecting  friction  losses,  any  injection  velocity 
becomes  theoretically  possible;  creating  a  possibility  of  per¬ 
forming  heat  loss  and  quality  calculations  for  an  injection 
velocity  that  is  impossible  due  to  excessive  friction  losses. 


' 
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The  mathematical  model  is  believed  to  be  reasonably 
comprehensive  and  any  difference  between  theoretical  and  actual 
results  could  presumably  be  traced  to  the  empirical  correlations. 
The  correlation  employed  to  determine  friction  losses  is  probably 
open  to  the  most  discussion ,  as  a  correlation  that  is  basically 
unreliable ( 10 )  is  subjected  to  extrapolations  and  interpola¬ 
tions.  A  straight  line  approximation  for  the  friction  factor 
versus  quality  relationship  simplified  the  calculations,  but 
presumably  led  to  errors  that  were  pronounced  for  qualities 
greater  than  0.95  or  less  than  0.4  (see  Figure  15  and  Figure 
17) .  Considering  the  correlation  and  the  approximations 
utilized,  it  would  be  optimistic  to  expect  a  completely  re¬ 
liable  temperature  profile.  But  until  complete  experimental 
or  field  data  are  available  to  evaluate  all  possible  two  phase 
flow  correlations,  it  is  felt  that  more  sophisticated  correla¬ 
tions  would  magnify  calculational  difficulties  without  necessar¬ 
ily  yielding  better  results.  Whereas  most  other  correlations 
(2,5,7,9,16)  were  derived  empirically  for  the  horizontal  flow 
of  oil  and  gas  mixtures,  the  relationship  employed  was  developed 
expressively  for  the  type  of  flow  encountered  in  two  phase  steam 
in je  ction . 

The  approximations  made  in  developing  the  f(t)  function 
(see  Appendix  I-C)  were  reasonable,  as  neither  the  change  in 
the  wellbore  heat  transfer  coefficient  nor  the  change  in  steam 
temperature  (see  Figure  3  and  Figure  11)  with  time  was  severe. 
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b )  Discussion  of  Results 

The  sample  cases  run  on  the  computer  yielded  results 
that  are  believed  to  be  reasonable  and  readily  explainable. 

A  simplified  computer  program,  neglecting  frictional 
effects  was  employed  in  an  attempt  to  duplicate  results  pub¬ 
lished  by  Satter(13)0  Figure  2  indicates  that  Satter  calcu¬ 
lated  heat  losses  that  were  approximately  four  percent  higher 
than  those  obtained  from  the  simplified  computer  program;  where¬ 
as,  Figure  3  indicates  a  difference  of  approximately  12°F  in 
the  calculated  casing  temperatures  «>  The  reason  for  the  dis¬ 
agreement  is  difficult  to  determine,  as  Satter  presented 
incomplete  data  by  neglecting  to  specify  values  for  emissivities 
and  annular  pressure.  He  also  neglected  to  specify  the  correla¬ 
tion  employed  to  determine  k  and  the  assumptions  made  in  cal¬ 
culating  f(t).  At  a  flowrate  of  3000  lb^/hr  the  effect  of 
friction  on  the  total  heat  loss  is  negligible,  as  shown  on 
Figure  2.  Increasing  the  flowrate  to  5000  lb  /hr,  the  neg- 
lection  of  friction  yields  heat  losses  that  are  five  percent 
higher.  The  frictional  effects  reduce  the  steam  temperature, 
subsequently  reducing  the  heat  losses  from  the  wellbore.  The 
effect  of  friction  on  the  total  heat  loss  is  not  pronounced , as 
temperature  changes  are  insignificant  at  shallow  depths.  If 
friction  is  neglected,  the  curvature  of  the  graphs  in  Figure  2 
is  solely  the  result  of  the  temperature  gradient  in  the  earth. 


' 


FIG.  2.  THE  EFFECT  OF  FRICTION  ON  THE  HEAT  LOSS 
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FiG.  3.  THE  EFFECT  OF  FRICTION  ON  THE  CASING  TEMPERATURE 
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Figure  4  represents  an  attempt  to  evaluate  the  model, 
specifically  the  two  phase  flow  relationship,  by  calculating  a 
temperature  profile  and  matching  it  to  an  observed  temperature 
profile.  All  the  known  data  (flowrate,  geothernal  gradient 
and  injection  pressure)  were  used,  while  reasonable  values 
were  assumed  for  the  rest  of  the  required  data.  The  computer 
results  for  two  different  injection  qualities  are  plotted  in 
Figure  4,  along  with  the  measured  field  data (3).  The  close 
agreement  obtained  for  an  injection  quality  of  0.9  does  not 
verify  the  reliability  of  the  model,  as  too  many  parameters, 
including  tubing  diameter  were  assumed.  However,  the  field 
data  presents  a  temperature  profile  that  is  similar  in  shape 
to  those  calculated,  which  would  indicate  that  the  model  is  a 
step  in  the  right  direction..  The  temperature  decrease  per 
foot  of  wellbore  increases  from  approximately  0.02  at  the  sur¬ 
face  to  approximately  0.04  at  the  sandface  (for  an  injection 
quality  of  0.9).  This  increase  is  due  to  lower  temperatures, 
higher  heats  of  evaporation  and  lower  heat  losses  combining 
to  increase  the  volume  of  vapor  flowing,  subsequently  increas¬ 
ing  the  friction  losses.  Increased  friction  losses  in  turn 
result  in  increased  temperature  decrements. 

The  effect  of  tubing  size  on  the  steam  temperature 
profile  is  shown  in  Figure  5.  Since  the  pressure  drop  for 

5 

single  phase  gas  flow  is  approximately  proportional  to  1/r^  , 
a  temperature  profile  that  is  very  sensitive  to  the  fractional 
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change  in  tubing  diameter  would  similarly  be  expected  for  two 
phase  flow.  Figure  5  verifies  this  hypothesis  as  multiplica¬ 
tion  of  the  tubing  diameter  by  1*125  (from  2"  to  2 V )  reduces 
the  temperature  drop  by  63  percent  (from  90°F  to  33°F)  ;  where¬ 
as^  multiplication  of  the  tubing  diameter  by  1.09  (from  2-3/4" 
to  3")  decreases  the  temperature  drop  by  35  percent  (from  10°F 
to  6^°F) . 

The  effect  of  tubing  size  on  the  quality  profile  is 
shown  in  Figure  6.  For  high  velocity  flow,  as  would  be  en¬ 
countered  in  small  diameter  tubing,  steam  temperature  changes 
would  be  high.  Thus,  with  increasing  depths  the  heat  losses 
will  decrease  and  the  heat  of  vaporization  will  increase, 
resulting  in  a  lower  fractional  condensation  or  quality  change. 
It  would  be  possible  to  have  the  quality  increase  with  depth 
if  the  loss  of  potential  energy  is  larger  than  the  heat  loss 
to  the  formation.  If  friction  were  neglected,  tubing  size 
would  not  have  any  effect  on  the  quality  profile.  Figure  6 
indicates  that  this  assumption  could  lead  to  considerable 
error  by  presenting  a  quality  change  of  0.12  for  two-inch 
diameter  tubing  and  a  quality  change  of  0.18  for  three-inch 
tubing . 

Figure  7  represents  a  cross-plot  of  sandface  pressure 
versus  tubing  radius  and  exemplifies  the  effect  of  tubing 
size  on  friction. 

The  total  heat  losses  are  lower  at  high  injection 
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rates,  as  indicated  in  Figure  9,  due  to  larger  steam  tempera¬ 
ture  changes  arising  out  of  increased  frictional  effects.  The 
sandface  quality  is  higher  for  high  injection  rates,  since  more 
available  vapor,  lower  heat  losses  and  larger  enthalpies  of 
vaporization  combine  to  cause  less  fractional  condensation. 

If  the  effect  of  friction  was  neglected,  the  curves  in  Figure  8 
would  be  expected  to  approach  an  asymptote  as  the  heat  loss 
per  pound  of  steam  flowing  becomes  negligible  compared  to  the 
potential  energy  loss. 

The  effect  of  injection  temperature  on  the  quality  pro¬ 
file  is  shown  in  Figure  10.  Since  the  heat  of  vaporization  and 
the  friction  losses  decrease  with  temperature  and  the  heat 
losses  increase  with  temperature,  higher  injection  temperatures 
would  be  expected  to  yield  increased  quality  changes.  Figure 
10  verifies  this  hypothesis. 

The  specific  volume  of  saturated  water  vapor  decreases 
with  temperature.  Thus,  higher  injection  temperature  would 
result  in  lower  velocities,  leading  to  lower  friction  losses. 
Figure  11  indicates  a  sandface  temperature  of  470°F  for  an 
injection  temperature  of  560  F  after  10,000  days.  For  an 
injection  temperature  of  640°F  the  sandface  temperature  be¬ 
comes  625°F.  The  bottomhole  temperature  decreases  with  time, 
as  the  volume  of  vapor  flowing  increases  with  time.  Figure  11 
also  indicates  that  the  change  in  steam  temperature  with  time 
is  relatively  small. 
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FIG.  11.  INJECTION  TEMPERATURE  VERSUS  SANDFACE  TEMPERATURE 
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FIG.  12.  THE  EFFECT  OF  INJECTION  TEMPERATURE  ON  THE  HEAT  LOSS 
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Figure  12  presents  the  effect  of  injection  temperature 
on  the  rate  of  heat  loss.  As  expected,  the  heat  loss  increases 
with  injection  temperature.  Figure  12  also  indicates  the 
effect  of  time  on  the  rate  of  heat  loss.  As  the  heat  losses 
decrease  by  27  percent  in  ten  days  (from  1  day  to  10  days  for 
an  injection  temperature  of  560°F) ,  it  is  apparent  that  an 
unsteady  state  solution  is  definitely  required.  All  the 
curves  should  and  appear  to  converge  to  zero  at  an  infinite 
time . 

The  effect  of  changing  various  parameters  on  the  rate 
of  heat  loss  is  shown  in  Figure  13,  It  is  interesting  to  note 
that  the  addition  of  a  second  casing  increases  the  heat  losses 
to  a  low  thermal  conductivity  formation.  The  second  casing 
was  added  outside  the  first  casing,  resulting  in  a  larger 
wellbore.  It  appears  that  the  resistance  of  the  annulus  is 
slightly  lower  than  the  resistance  of  the  formation  it  replaces. 
For  high  thermal  conductivity  formations,  the  addition  of  the 
second  casing  decreases  the  rate  of  heat  loss.  Considering  a 
formation  with  a  thermal  conductivity  of  1.4  Btu/hr . ft . °F,  the 
addition  of  a  10%"  casing  reduced  the  rate  of  heat  loss  after 
10  days  by  20  percent  (see  Figure  13) .  The  reduction  is  even 
larger  at  shorter  times. 

From  Figure  13  it  becomes  evident  that  radiation  in 
the  annulus  is  a  very  important  mode  of  heat  transfer,  since 
decreasing  the  tubing  emissivity  from  0.8  to  0,2  reduced  the 
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FIG.  13.  THE  EFFECT  OF  PARAMETER  VARIATIONS  ON  HEAT  LOSS 
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rate  of  heat  loss  by  30  percent  (after  10  days) .  Reduction  of 
the  earth  thermal  conductivity  from  1.4  Btu/hr.ft.°F  to 
0.4  Btu/hr.ft.°F  decreased  the  rate  of  heat  loss  by  47  percent 
(after  10  days).  Thus,  the  earth  thermal  conductivity  and  the 
emissivities  become  important  parameters  in  determining  a 
representative  rate  of  heat  loss  from  the  wellbore. 

The  discussion  of  the  results  may  be  summarized. 

1)  Higher  velocities  resulted  in  higher  friction 

losses . 

2)  Higher  friction  losses  resulted  in  lower  steam 
temperatures . 

3)  Lower  steam  temperatures  resulted  in  lower  heat 
losses,  higher  heats  of  evaporation  and  higher  friction  losses, 

4)  Lower  heat  losses  resulted  in  smaller  quality 

changes . 

5)  Higher  heats  of  evaporation  resulted  in  smaller 
quality  changes. 

6)  Smaller  quality  changes  resulted  in  increased 
friction  losses. 

Due  to  the  large,  number  of  variables  involved,  it  was 
difficult  to  evaluate  the  effect  of  every  variable  on  the 
results.  Thus,  the  graphs  presented  in  this  discussion  are 
meant  to  be  illustrative  only,  and  to  apply  the  proposed  model, 
the  computer  program  would  have  to  be  run  using  the  specific 
data  applicable  to  the  problem  to  be  solved. 
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CONCLUSIONS 

The  following  conclusions  were  made  on  the  basis  of 
the  results  of  the  sample  computer  runs. 

1)  Consideration  of  friction  led  to  lower  heat  losses, 
but  the  effect  was  relatively  minor;  whereas,  the  quality  and 
the  temperature  profiles  were  greatly  affected  by  friction. 

2)  The  rate  of  heat  loss  may  be  reduced  by  the 
addition  of  a  second  casing  or  the  reduction  of  the  emissivities . 
The  effectiveness  of  the  addition  of  a  second  casing  larger 
than  the  first  casing,  is  governed  by  the  thermal  conductivity 
of  the  earth.  If  the  thermal  conductivity  is  low,  the  second 
casing  may  increase  the  rate  of  heat  loss. 

3)  The  earth  thermal  conductivity  is  a  very  important 
parameter  in  the  determination  of  the  rate  of  heat  loss. 

4)  The  shapes  of  the  calculated  temperature  profiles 
were  comparable  to  an  observed  temperature  profile. 

5)  The  effect  of  the  time  on  the  steam  temperature 
and  the  overall  wellbore  heat  transfer  coefficient  was 


relatively  minor. 
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RECOMMENDATIONS 

The  assumption  of  a  constant  wellbore  steam  temperature 
is  recommended  for  cases  where  only  the  rate  of  heat  loss  is 
desired.  To  ensure  that  friction  losses  are  not  excessive, 
a  rough  check  on  the  possible  friction  losses  would  be  advis¬ 
able.  The  assumption  would  simplify  the  calculations  con¬ 
siderably,  and  the  error  introduced  would  act  as  an  inherent 
safety  factor.  To  calculate  representative  heat  losses,  care 
must  be  exercised  in  arriving  at  a  value  for  the  thermal  con¬ 
ductivity  of  the  earth.  If  conductivity  data  is  available  as 
a  function  of  depth,  and  considerable  variations  are  indicated, 
it  might  be  advantageous  to  divide  the  formation  into  strata 
of  constant  physical  properties.  For  situations  involving 
small  steam  and  earth  temperature  changes  with  depth,  and  the 
total  heat  loss  is  the  desired  result,  an  average  value  for 
the  thermal  conductivity  is  believed  to  yield  sufficient 
accuracy . 

To  calculate  reasonably  reliable  values  for  the  sand- 
face  pressure  and  quality,  the  effect  of  friction  must  be  con¬ 
sidered.  Due  to  rapid  error  propagation  resulting  from  the 
relationship  between  pressure  and  friction  losses,  a  very  good 
two  phase  flow  correlation  is  required  to  obtain  reliable  sand- 
face  pressures.  Until  a  correlation  is  derived  that  has  been 
substantiated  under  conditions  similar  to  field  conditions,  it 
is  recommended  to  employ  caution  in  using  calculated  sandface 
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temperatures,  as  the  error  could  conceivably  be  considerable. 

The  computer  program  as  it  is  included  in  Appendix  II, 
could  readily  be  modified  to  account  for  more  resistances  in 
the  wellbore  heat  transfer  coefficient  and  more  elaborate 
relationships  between  air  properties  and  temperature.  These 
changes  are  not  recommended,  as  it  is  believed  that  the  small 
gain  in  accuracy  will  not  justify  the  increased  computational 
time.  If  at  any  time  the  two  phase  flow  correlation  is  sub¬ 
stantiated,  it  may  be  advantageous  to  employ  a  higher  order 
polynomial  for  the  relationship  between  friction  factor  and 
quality. 

Of  the  means  of  decreasing  the  heat  losses  from  the 
wellbore,  the  reduction  of  the  emissivities  is  recommended 
above  the  addition  of  a  second  casing.  Whereas, the  reduction 
of  the  emissivities  always  reduces  the  heat  losses,  the  addi¬ 
tion  of  a  second  casing  (larger  than  the  first  casing)  may  in¬ 
crease  the  heat  losses.  The  extent  and  under  which  circum¬ 
stances  the  emissivities  are  to  be  reduced  becomes  a  matter  of 


economics . 
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NOMENCLATURE 

=  heat  capacity  of  the  annulus  air  (Btu/lbm°F) 

=  dH/dT  of  water  at  its  bubble  point  (Btu/lbm°F) 

=  two  phase  friction  factor  as  correlated  in 
Figure  14 

=  time  function  as  defined  in  the  Literature  Review 
=  Grashof  number 
=  enthalpy  (Btu/lbm) 

=  radiation  heat  transfer  coefficient  as  defined  by 
equation  B2 

2 

=  film  coefficient  on  inside  of  tubing  (Btu/hr.ft.  °F) 

=  enthalpy  of  vaporization  (Btu/lb  ) 

=  integral  as  defined  by  equation  A19 

=  integral  as  defined  by  equation  C13 

=  mechanical  equivalent  of  heat  (778  ft-lb^/Btu) 

=  thermal  conductivity  of  the  annulus  air  (Btu/hr . ft . °F) 

=  effective  thermal  conductivity  of  the  annular  air 
space  (Btu/hr . ft . °F) 

=  thermal  conductivity  of  the  earth  (Btu/hr . ft . °F) 

=  thermal  conductivity  of  the  tubing  (Btu/hr . ft . °F) 

=  thermal  conductivity  of  the  casing  (Btu/hr . ft „ °F) 

=  thermal  conductivity  of  the  second  casing 

(Btu/hr. ft . °F ) 

=  molecular  weight  of  air  (28.9  lb  /lb. mole) 

=  steam  pressure  (psia) 

=  air  pressure  in  the  annulus  (psia) 

=  Prandtl  number 
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average  rate  of  heat  loss  (Btu/ft.hr.)  as  defined 
by  equation  A4 

heat  gain  of  the  steam  (Btu/lb  ) 

m 

Gas-law  Constant  (10.73  lb^ft "Vin^lb .mole  °R) 

radius  of  the  wellbore  (ftO 

vapor  phase  Reynolds  number 

inside  radius  of  the  tubing  (ft) 

outside  radius  of  the  tubing  (ft) 

inside  radius  of  the  casing  (ft) 

outside  radius  of  the  casing  (ft) 

inside  radius  of  the  second  casing  (ft) 

outside  radius  of  the  second  casing  (ft) 

time  (hrs . ) 

temperature  of  the  outside  of  the  last  casing  (°R) 

dimensionless  time  as  defined  by  equation  C9 

earth  temperature  (°R) 

steam  temperature  (°R) 

steam  temperature  (°K) 

temperature  (°R)  at  r^ 

velocity  (ft/sec)  or  dummy  variable 

overall  heat  transfer  coefficient  based  on  r 
(Btu/hr.ft.  °F ) 

dimensionless  coefficient  as  defined  by  equation  C8 

3 

specific  volume  (ft  /lb  ) 

steam  flowrate  (lb  /hr.) 

m 

quality  (fraction) 
depth  (ft) 


a 


Y 


A 


t 


A 

c 


A 


P 


a 


u  G 
£1 
e  2 


e  3 


E  4 


2 

=  thermal  diffusivity  of  the  earth  (ft  /hr. ) 

=  parameter  as  defined  by  equation  A10 

=  temperature  drop  across  the  tubing  (°F) 

=  temperature  drop  across  the  casing  (°F) 

=  temperature  drop  across  the  second  casing  (  F) 

=  density  (lb  /ft3) 

m 

—  8 

=  Stefan-Boltzman  constant  (0.1713  x  10  Btu/hr.ft 

=  viscosity  of  the  annulus  air  (lb  /hr. ft*) 

m 

=  viscosity  of  saturated  water  vapor  (lb  /hreft0) 

=  emissivity  of  the  outside  of  the  tub  ng 
-  emissivity  of  the  inside  of  the  casing 
=  emissivity  of  the  outside  of  the  casing 
=  emissivity  of  the  inside  of  the  second  casing 


Note  : 

Any  of  the  above  symbols  subscripted  with: 

1)  a  G,  represents  the  vapor  phase,  e.g.  pG 

2)  an  L,  represents  the  liquid  phase,  e.g. 

3)  a  number,  represents  the  point  in  the  tubing,  e.g.  y^ , 

or  the  annular  space  number,  e.g.  k  ,  depending  on 

C2 

context. 
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DERIVATION  OF  EQUATIONS 
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I -A.  THE  FLUID 

(Calculation  of  Quality  and  Depth) 

Writing  a  differential  energy  balance (14)  for  the  steam 
flowing  down  the  tubing  yields: 

g  dz  udu  dW 


dH 


g  J 


J9, 


=  dq 


w 


(Al) 


As  there  is  no  shaft  work,  dW  =0. 

'  s 

Defining  q^  as  the  rate  of  heat  loss  per  foot  of  pipe  yields 


dq 


w 


dz 


W 


st 


Thus,  equation  Al  reduces  to: 


dH  + 


udu 


Jg. 


=  (■ 


g  1 


-)  dz 


(A2) 


W 


st 


Assuming  a  depth  interval  small  enough  to  assume  linear  casing 
and  earth  temperature  profiles,  and  integrating  equation  A2 
from  point  1  to  point  2  yields: 

i  g  i  q„ 


H2  -  H1  + 


2  Jg 


,  2  2 . 
(u2  -  Ul  ) 


( 


g  J 
Jc 


W 


*>  (Z2  ~  zq) 


St 


(A3) 


f(t)  is  defined  to  make  the  following  equation  correct (11) 


q  dz 
nr 


2*  kE  (TC  -  V 

f  (t) 


dz 
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Thus  if 


and 


=  a  +  bz 

£j 


T  =  0  -f  fz 


2tt  k. 


(e  +  fz  -  a  -  bz)  dz 


z2  -  zi  zi 


Defining : 


T  =  e  +  fz, 
ci  1 


T  =  a  +  bz.. 

Ei  1 


T  =  e  +  fz0 
c2  2 


'E. 


a  +  bz. 


and  assuming  f(t)  is  constant  over  the  interval  z,  to  z„  yields 


it  k. 


f  (t) 


(T  -  T  )  +  (T  +  T  ) 

-  C1  L1  C2  b2  - 


(  A4 ) 


If  the  vapor  phase  may  be  assumed  to  be  the  principal  factor  in 
determining  kinetic  energy  effects,  obtain: 


W, 


y  W  ,  =  -rr 

J  st 


rl  p g  UG  3600 


and 


W 


st 


u  =  u. 


Trr-^  *  3600 


y  vr 


The  evaluation  of  the  kinetic  energy  contribution  leads  to 


,2  2,  e  ,  2  2  2, 
- -  (u„  -  u,  )  =  f  (y  v  -  y,  v  ) 

2g  J  J2  G, 


(A5 ) 
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where 


W 


st 


2gcJ 


irr^  3600 


Defining  H  as  a  reference  enthalpy,  can  write 


II  „  =  II  + 


T  T 

s,  : 

f  C  ,  dT  +  f 
ip  sL  s  ^ 


C  _  dT 
sL  s 


+  H  y  0 
v2  2 


'1 


Thus  : 


H 


1 


H  + 
o 


T 


T 

o 


+ 


y 


2 


y 


l 


For  reasonable  temperature  changes,  C  _  is  nearly  constant  and 

S  i-j 

thus  obtain: 


H 


2 


m 


)  +  y 


o  H 
2  v 


2 


(A6) 


Substituting  equations  A5  and  A6  into  equation  A3  and  solving 
for  y£  yields: 


{yx(Hv  +  y^v,,  ) 


1  1  G- 


C  _  (T  -  T  )  +  (- - 

sL  s~  s,  Jg 

2  1  ^  c 


q. 


w 


st 


)  ( z 2  -  Zp) 


(H 


v. 


+  y2  fl  VG„  > 


Equation  A7  is  implicit  in  y2 .  Since  the  kinetic  energy  con¬ 
tribution  is  generally  quite  small  and  a  trial  and  error  solu¬ 
tion  is  required  for  the  overall  solution,  the  y9  on  the  right 


} 

-  (A7) 


--  -  —  / 

l. 


s" 


...  - 


X sm  2  sJ’ix/p  £ll£n  >t  «>p 
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hand  side  will  have  little  or  no  effect  on  the  rate  of  conver¬ 
gence;  whereas,  an  explicit  equation  would  require  the  solution 
of  a  quadratic,  which  in  turn  would  be  time  consuming. 

The  equation  describing  the  pressure  profile  for  the 
vertical  flow  of  condensing  steam  down  a  circular  pipe,  is 
taken  from  Perry (10)  p.  5  -  42. 

It  may  be  written  as: 
dP 
dz 


4fG  PG  UG 


(A8) 


where 


2g  2rn 
Jc  1 


fG  =  F(REG’Y) 


(A9) 


(The  graphical  relationship  for  f^  was  taken  from  Perry (10) 
p.  5-40,  enlarged,  and  included  in  Appendix  III-A  as 
Figure  14 . ) 

(1  -  y)  W 


Y 


st 


2"  Pl  rl 


(A10) 


u. 


y  W 


st 


p  it  r  n  3600 


G  1 


R 


EG 


2pG  ug  2y  Wsfc 


■nr,  p 


1  MG 


Equation  A8  may  be  rewritten  as: 
dP 


v. 


=  Xf  y  dz 


(All) 
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whe  re 


A  =  (. 


W 


st 


-1 


12  *  36 00  *  tt 


■) 


g  rn 

i 


If  the  temperature  change  in  the  steam  is  not  excessive,  average 
values  may  be  assumed  for  u_  and  p_ .  Thus,  for  every  run,  a 
plot  of  quality  versus  two  phase  friction  factor  may  be  obtained, 


as  : 


where 


REG  =  YaY 


=  (X  -  y) • Y 


b 


2W 


st 


Trr,  u  „ 
I  ^ 


(A12) 


W 


st 


2"rlPL 


(A13) 


It  was  found  that  for  relatively  high  Y  and  Y,  ,  as 

a  d 

would  be  encountered  in  a  field  case;  a  straight  line  yielded 
a  reasonable  fit  for  f_  vs  y  in  the  interval  0.3  <y  <0.95 

YJ 

(see  Figure  15  and  Figure  17) .  The  coordinates  of  this  line 

were  nearly  independent  of  Y  and  Y,  (see  Figure  17).  Thus: 

a  d 

f„  =  a+  b.y  (A14 ) 

G  y  y  y 

From  Keenan  and  Keyes (8)  p.  14  and  p.  15  obtain: 


-  *<V  ' 


P 


srlJ 


i  j 3^'  {LfOJDfiB  3c  luov; 
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Therefore , 


dP  =  (T  )  dT 

s  s 


(A15) 


vr  =  iKT  ,P)  =  ip(  T  ,  $  (T  )  )  (A16 ) 

u  s  s  s 

Substituting  equations  A14,  A15  and  A16  into  equation  A8  , 
and  integrating  from  point  (1)  to  point  (2)  yields: 


T 

I 


T 


♦'  (Ts> 

+  (Tg,  *(Ts)) 


(a  +  b  y) 

y  y 


2  , 

y  dz 


If  the  change  in  quality  over  the  interval  to  z ^  is  not 
excessive,  use  of  an  average  quality  over  the  interval  is 
believed  to  result  in  sufficient  accuracy. 

Thus  : 


where 


-v  -2 

y)  y 


I 


T 


4>'  <v 

- - -  dT 

<MT  ,  <KT  )  ) 


(A18 ) 


(A19) 


y 


yi  +  y- 


Equation  A7  may  be  substituted  into  equation  A18,  resulting  in 

one  equation  with  two  variables  (T  and  z?).  Specification 

s2  ^ 

of  T  allows  the  determination  of  z?  by  a  trial  and  error 
s2 


(A17) 


' 
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solution  of  the  implicit  equation.  Newton's  method  was  applied 
for  faster  convergence. 


Fz  ■  Az 


—2  — 
X  y  ( a  +  b  y ) 

Y  y  1 


=  0 


(A20) 


Az(n+1)  =  Az (n) 


(n) 


(^Z)  (n) 


dAz 


(A21) 


dF  I  (2a  +  3b  y) 

_ £  =  i  +  T  y  y 

dAz  xy  (a 


+  b  y) 

y  y 


dy 

dAZ 


(A22) 


Equation  A7  was  used  to  obtain: 


dy 


dAz 


1  dy. 


g  J 
Jc 


W 


st 


2  dAz 


2  (H  +  y0  fn  ) 


v. 


2  n  G. 


(A23) 


The  substitution  of  equations  A22  and  A23  into  equation  A21 
provides  a  procedure  to  obtain  a  Az  to  satisfy  equation  A20. 

The  equation  for  the  vapor  pressure  of  steam  is  given 
in  Keenan  and  Keyes  (8)  as: 

-2 . 30  3u^, 


P  =  4>(Ts) 


P  e 
c 


4>'  (T  ) 

<3 


2.303  c*>(Ts)  u^ 


where 


u. 


sk 


3  4 

a  +  bx  +  Cx  +  ex 
P  P  P_ P 

(1  +  d  x) 

P 


•  •  '  • 


'-do  ~'1  u-j.  >bw  V  •  c  .  .  pi 


2  lO  .-r.it..'  ■>  -v'.  n«:  is  ::  r  ia;.  .  '? 
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du. 


du. 


t 


E 


E 


u. 


E 


DT 


sk 


dx 


3 


+  5e  x 
P 


4 


(T  -  x)  (1  +  d  x) 
c  p 


x 


T 


647. 27°K 


c 


a  =  3.3463130 

P 

b  =  4.14113  x  10"2 

P 

C  =  7.515484  x  10~9 

P 

d  =  1.3794481  x  10"2 

P 

e  =  6.56444  x  10"11 

P 

P  =  critical  pressure  of  the  steam  (psia) 

The  determination  of  v_  is  straight  forward.  Equation  13 

page  15  in  Keenan  and  Keyes (8),  was  applied  with  the  use  of  a 

1  3 

conversion  factor  to  change  the  units  to  ft  /lbm. 


IT  was  subsequently  evaluated  numerically  using  a  nine 


point  Gauss-Quadrature  integration  formula. 


t  IT  was  evaluated  in  degrees  Kelvin,  since  the  equation  for 
vr  and  P  were  given  in  degrees  Kelvin. 


ji  :iq  J  fi' 

.  B  -  um  ■  v' 


-r>r 

P-  1 
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I-B .  THE  WELLBORE 

(Calculation  of  Casing  Temperature  and  UQ) 

1)  Single  Casing 

At  steady  state,  with  no  heat  source  or  sink,  the  heat 
transferred  through  every  resistance  in  the  wellbore  must  be 
equal (6,  12).  Thus  referring  to  Figure  1  obtain: 


<$q  =  21Tr.dz  h  (T  T.  ) 

j.  s  s  J. 


=  2vdz  k. 


T  -  T 
1  2 


In 


=  2Trr2dzf12  (T2  +  T3)  (T22  +  T32)  (T2  -  T3) 


+  2ir  dz  k 


<T2  -  T3) 


In 


=  2frdz  k. 


T3  -  T4 


In 


Solving  for  the  temperature  differences,  adding  them 
and  subsequently  solving  for  6q,  yields: 


2 tt  dz  (T 


6q  = 


2irr/,dz-U  (T  -  T.)  = 
4  o  s  4 


T4 ) 


(BJ 


In 


In 


r,  h 
1  s 


r0h  + 
2  r 


In 


where 


f  12  cr  (T2  +  T3)  {T22  +  T32) 


(B2 ) 


12 


r~  1 

+  —  -  1) 


r3  S  2 


Solving  equation  B1  for  Uo  obtain: 


1 

U 


r.  In  — 
4  r. 


r  *  In  — 

4  r. 


rlhs 


r~h 
2  r 


Letting : 


4 


r ,  In  — 
4  r. 


r,  h 
1  s 


—  f  t  0a 

r4  12 


r  „  In  — 1 
4  r. 


Equation  B3  may  be  rewritten  as: 


In 


r ,  In  — 
4  r. 


(B3) 


U 


1 

c. 


r 


i  - 


C3C4 (T2  +  T3) (T22  +  T32)  +  C3C5kc  +  1 


CB4) 
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Where  kc  is  defined  as  the  effective  thermal  conductivity  of 
the  air  in  the  annular  space.  It  combines  natural  convection 
and  conduction.  From  Groeber(6)  page  317  obtain: 


k  =  0.11k  (Pr’Gr)0'29  (2000+  < Pr*Gr <  106) 

c  a 

k  =  k  (2000  >  Pr-Gr) 

a. 


(B5 ) 


Pr*  Gr  = 


7500(12)  Cp  (T2 
a 


T3)  Pa2M2(r3 


3  2 

y  k  tm  R 

a  a  M 


C6  CP  (T2  -  T3> 
a  _ 

y  k  T  3 
a  a  M 


(B6) 


where 


M 


T2  +  T3 


7500  (12)  3  Pa2  M2  (r3  -  r2)3  g 


R 


u  =  F (T  ) 
a  v 


S  -  F(V 

a 


all  functions  are  assumed  linear 


ka  =  F(Tm) 


Since  the  heat  transferred  through  the  wellbore  must  equal 
the  heat  entering  the  formation  obtain: 


t  the  relationship  is  extrapolated  from  the  specified  minimum 
of  6000  to  2000.  The  error  introduced  is  negligible. 


f 
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2"k 


2*r4  (Ts  '  V  Uodz  = 


E 


f  (t) 


(Tc  -  Te)  dz 


Solving  for  T  yields: 

c 


=  T 


kE  (Ts  -  V 

kE  +  Uor4  f(t> 


or 


T  -  T  + 
c  s 


kE  (Ts  -  TE) 

kE  +  Uor4  f(t> 


=  0 


(B7 ) 


Assuming  the  temperature  drops  across  the  tubing  and  the 
casing  are  fai  y  constant,  may  write: 


T3  =  To  +  Ac 

Thus,  Tc  remains  the  only  variable  in  equation  B4. 

Substitution  of  equation  B4  into  equation  B7  results  in  one 

equation  in  T  .  Due  to  the  complicated  nature  of  the  equation, 

o 

a  trial  and  error  solution  appears  to  be  necessary.  For  a 
fast  and  reasonably  reliable  convergence  procedure,  Newton's 
method  was  employed. 

Thus  : 


T(n+1)  =  T(n) 

c  c 


(n) 


(  T  ) (n) 


dT 


. 
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dF 


T 


dT 


m  kE  (Ts  -  V  r4  f(t> 


dT. 


(kE  +  °o  r4  f  (t)) 


dT. 


The  differentiation  of  equation  B3  yields: 
dU 


dT. 


C.  (T  2  +  2T0Tq  +  3T  2 )  +  C-k ' 

4  A  A  i  D  DC 

Cc3*C4(T2  +  T3)(T22  +  T32)  +  C3*C5kc  +  fl2 


where 


k'  = 
c 


dk 


dT. 


k' 

a 


-  a 


0.2  9  d(Pr-Gr) 


Pr  *Gr 


dT. 


(B8 ) 


d  (Pr*Gr) 


dT. 


C ' 


k' 

a 


»  — 


C6 

^  2Cp 

a 

2k  v 
a 

T  3 
a  M 

-(T2 

-  T3> 

r3 

cp 

aLTM 

dCP 

a 

i  dcp 

a 

dT3 

2  dTA, 

M 

dk 

a 

1 

dk 

a 

dT3 

2 

dTM 

1 

dya 

dT3 

2 

dTM 

T3)  2CP 


2k' 

a 


u 


(B9) 
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After  Tc  has  been  calculated  by  the  prescribed  or  other 
iterative  procedure,  A  and  A  may  be  calculated  and  T 

"C  C  G* 

redetermined.  Comparison  of  the  two  values  so  obtained  will 
indicate  whether  or  not  another  iteration  is  required. 


A 


t 


A 


c 


kE(Tc  '  V 

f  (t) 


r,  h 
1  s 


k  (T 
E  c 


V 


In 


f  (t) 


2 )  Double  Casing 

The  treatment  of  the  double  casing  situation  is  similar 
to  the  single  casing  case.  Thus,  emphasis  will  be  placed  on 
the  differences  in  the  two  problems.  Following  the  same  pro¬ 
cedure  as  in  Part  (1)  obtain: 


1 

U 

o 


+ 


k 


+ 


h 


+ 


c2 


rr  In  — 
6  r 


(BIO) 
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1 


U 


oM 


r,  h 
1  s 


+ 


r ,  In  ~— 
4  r. 


1 


+ 


— -  h  + 

r4  r 


r .  In  — 
4  r. 


r  ,  In  — 
4  r. 


(Bll) 


where  UqM  is  defined  as  the  overall  heat  transfer  coefficient 
between  the  first  casing  and  the  steam,  based  on  the  outside 
area  of  the  first  casing. 

Equation  Bll  may  be  rewritten  employing  the  definitions 
of  Part  ( 1)  . 


U 


oM 


1 

C. 


1 


1  - 


C3  *  C4  <T2  +  V  <T22  +  T32)  +  C3'C5'kc  +  1 


Similarly 


U 


o 


C 


10 


FI  +  F2 


C x Q  9  F 1  *  F 2  +  FI  +  F2 


(B13) 


where 


r6 

.JL 

r  in  — 
6  r^ 

r6 

In 

r  3 

+ 

r r  In 
0 

^10 

o 

i — i 

u 

1 

kl 

k2 

k3 

FI 

r — 1 
rH 

u 

II 

(T2 

+  T3>(T22 

+  T32> 

+  C12 

k 

c 

F2 

c 

^13 

<T4 

+  T5)(P2 

+  t52) 

+  C1 4 

k 

C  9 

C 


11 


— —  f ,  ~  a 
r6  12 


(B12) 
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k  may  be  evaluated  by  equation  B5,  while  k  may  be  calculated 
c  c2 

similarly  using  the  properties  of  the  second  annulus 
T4  +  T5 

(T  =  - - -  ,  etc.)  . 

2  2 


At  steady  state,  the  heat  conducted  through  any  part 
of  the  wellbore  must  be  equal  to  the  heat  absorbed  by  the 
formation.  Thus: 


2lTkE(Tc  ~  V  dz 

f  ( t) 


2vrr  U  (T  -  T  )  dz 
6  o  s  c 


(B14 ) 


2irk_,  (T  -  Tr, )  dz 
E  C  E 

f  (t) 


2irr.U  .. 
4  oM 


~  )  dz 


(B15) 


Equation  B15  may  be  rearranged  to  yield: 


+ 


(B16 ) 
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T 


c 


r4  f(t) 


The  solution  of  equation  B14  for  T  yields 


T  =  T  + 
c  s 


kE  ^TE 


T  ) 
s 


kE  +  rg  f(t>  Uo 


or 

kE  (Ts  -  V 

F  ?  =  T  -  T  +  - - - - - - -  =  0  (B17) 

kE  +  r6  f <t)  Uo 

It  is  now  possible  to  substitute  equations  B12,  B13  and  B16 
into  equation  B17  and  derive  one  implicit  equation  in  (or  T^) . 
Newton's  method  was  used  as  a  basis  for  convergence. 


m  (n+1) 

3 

m  (n) 

3 

F(n) 

T2 

(dFT2} (n) 
dT3 

(B18) 

Assuming 

T2 

=  T 

s 

At 

T4 

=  T  - 

3 

A 

c 

T  c. 

=  T  + 

c 

A 

C2 

dFT2 

dT 

c 

kE  <Ts 

-  Te)  r6  f(t) 

dU 

o 

dT3 

dT3 

(kE  + 

r  r  f  (t  )  U  )  ^ 

6  o 

dT3 
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The  differentiation  of  equation  B16  yields: 


dT 


dT. 


r4 


"E 


dU 


-  U  +  (T  -  TJ 
oM  s  4 


OM 


dT. 


The  differentiation  of  equation  B12  yields 
dU 


oM 


C4(T22  +  2T2  T3  +  3T32)  +  C5  k; 


dT3  Qc3-c4  (t2  +  t3)(t2  +  t3  )  +  c3-c5kc  +  ll 

where  k'  is  defined  by  equation  B8. 
c 

The  differentiation  of  equation  BX3  yields: 


dU 


o 


dT. 


10 


(FI  +  F2)  [_  C1Q‘  ( FI  ’  F2  '  +  F2-F1') 


L  C1Q  *F1 -F2  +  FI  +  F2  J 


T2" 


FI'  +  F2' 


(C10*F1*F2  +  FI  +  F2) 


where 


FI'  = 


dFl 


dT. 


=  C11(T22  +  2T2  T3  +  3T32)  +  C12ki 


F2 


dF2  ~  7 

—  =  C13(T5  +  2T4  T5  +  3T4  >  +  C14k;, 

dT~  2 


+  [c13(T42  +  2T4T5  +  3T52)  +  C14k^  J 


dT 


2  dT. 


3k 


k' 

c. 


3T . 


(T^  is  constant) 


3k 


k' 

c. 


3T, 


(T^  is  constant) 


+  Fl'  +  F2'J 


After  a  value  for  Tc  has  been  obtained,  better  values 

for  A  A  and  A  may  be  calculated. 

c2 


A  =  t  -  T 
t  s  2 


A  =  T  _  -  T  . 
c  3  4 


=  TV  -  T 
c2  5  c 


k—,  (T  - 
E  c 

V 

1 

( 

f  (t ) 

\ 

r,  h 

1  s 

r  4 

*■* 

M 

HJ 

0 

1 

te) 

In  — 
r  3 

f  (t) 

k2 

kE^Tc 

te> 

t  r6 
In  — - 

r5 

f  (t) 

k3 

In 


k. 


It  was  found  that  A  ,  A 

t  '  c 

did  not  affect  the  value  for  T  . 

c 


and  A  were  quite  small,  and 
c2 

Therefore,  no  iterations 


A  and  A 
c  c. 


were  necessary  to  correct  the  values  for  A^, 


' 
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I-C.  THE  FORMATION 
(Calculation  of  f(t)) 


The  partial  differential  equation  describing  conduction 
in  cylindrical  coordinate  for  a  constant  thermal  conductivity 
is (12)  : 


1  32T  1  3  3T  32T  1  3  T 

—  -  +  —  —  (r  — )  +  -  =  -  (Cl) 

r2  3  G2  r  3r  3r  3z2  a  3t 

Assuming:  radial  symmetry 

no  axial  (vertical)  conduction 

Equation  Cl  reduces  to: 

3  2T  1  3T  1  3T 

- + - =  -  —  (C2 ) 

3r2  r3r  a3t 


Defining : 

T*  =  T  -  T  Therefore,  T  =  T*  +  T 

s  '  s 

T  *  =  T_,  -  T 

O  E  s 

The  substitution  of  the  expression  for  T  into  equation  C2 
yields : 

32T*  1  3T*  1  3T*  3T 

- + - =  -  ( -  +  — -)  (C3 ) 

3r2  r  3r  a  3t  3t 


3T 

Since  -  is  generally  quite  small, 

O  t 

zero,  whereby  equation  C3  reduces  to 

32T*  1  3T*  1  3T  * 

3r2  r  3r  a  3t 


it  was  assumed  to  be 


(C4) 


' 
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Equation  C4  is  solved  in  Carslaw  and  Jaeger (4)  page  337 
employing  Laplace  transforms  and  a  contour  integral.  The 
solution  is  subject  to  the  following  initial  and  boundary  con¬ 
ditions  : 

Initial  condition: 

at  t  =  0  and  r  ^r  ,  T*  =  T  *  (i.e.  T  -  T_) 

c'  o  E 

Boundary  conditions: 

1)  Temperature  (T)  is  finite  at  an  infinite  radius 

2)  Heat  conducted  through  the  wellbore  is  equal  to 
the  heat  entering  the  formation.  In  mathematical  form  obtain: 


as 


2tt r  U  (T  -  T  )  dz  =  -  k„2ur 

cos  c  E  c 


3T 


3  r 


dz 


r=r 


or 


U  (T  -  T  )  =  k„ 
o  c  s  E 


3T 


3r 


r=r 


Thus 


3T ' 


3  r 


U 


rp  * 


r=r 


"E 


r^r 


(C5 ) 


3T*  3T 


and 


3r  3  r 


T  -  T  =  T* 
c  s 


r=r 


The  assumption  of  a  constant  heat  transfer  coefficient  and  the 
stated  boundary  and  initial  conditions  lead  to  the  following 
solution  of  equation  C4 . 


.  t 
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T*(r,t)  = 


U 

■2  r~  T  *  0 

k  o  00  2, 

E  r  -au  t 

- -  /  e 


7T 


U  U 

J  (ur)  (uY,(ur  )  +  Y  (ur  )  }  -  Y  (ur)  {uJ,  (ur  )  +  ~  j  (ur  )  } 

1  C  KE  °  C  °  1  c  kE  °  

{uJl(urc>  +  kT  Jo (urc' )2  +  {uYl(urc>  +  kT  Yo(urc) ,2 

iii  t 


du 


u 


(C6 ) 


where  Jq (ur) ,  (ur) ,  YQ(ur)  and  (ur)  are  all  Bessel  functions 
Since  the  heat  flux  at  the  wellbore-formation  interface  is  the 
desired  property,  may  write: 


Thus : 


3T 

3T* 

U 

_  o 

3r 

3r 

r=r 

c 

II 

n 

o 

X 

M 

3T 

u  9 

■2  (jT1 2)2  T  * 
ke  ° 

3  r 

r=r 

c 

TT 

T* 


r=r 


(from  equation  C5) 


/ 

o 


2. 

-otu  t 


Jo  (urc}  -  Yo(urc>  ujl(urc> 

{uJl(urc>  +  £  JQ(urc)}2  +  {uYl(urc)  +  ^  Vurc)}2 

iii  t 


du 

u 


But 


-Jo(urc)  YI(urc)  +  Jl(urc>  = 


tt  (ur  ) 
c 


* 


V  -  >.  -  i  <■  -  , 
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Thus  : 


3T 

4U  2t*  «> 

T  o  c 

i — ■ 

N 

3 

Q 

4-) 

1 

<D 

! 

j 

L 

3  r 

2  ^ 
it  r  o 

r=r  c 

c 

(o)J-  (a>)  +  UmJ  (a>)  }  2+  {  wY ,  ( a) ) 
1  To  1 

+  UTYo(“),2_ 

where  : 


do) 

OJ 

(C7 ) 


oj  =  ur 


U„ 


U  r 
o  c 


(C8 ) 


at 


'D 


(C9 ) 


From  the  definition  of  f(t)  obtain 

3T 


2irr  dz 
c  E 


3r 


r=r 


2*kE  (Tc  -  V 

f  (t) 


dz 


Solving  for  f(t)  yields 

f  (t)  = 


T 

E 


3T 
c  3r 


r=r 


From  equation  C5  obtain 


T  =  T  + 
c  s 


r  3T 
c 

UT  3  r 


Thus 


f  (t) 


T  -  T 
E  s 


3T 
c  3r 


r=r 

< 

1 

U 


T 


r=r 


(CIO) 


i 
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The  substitution  of  equation  CIO  into  equation  C7  yields 


IT 


f  (t) 


( - ) 

2Um 


1 

IL 


(Cll) 


where  : 


h  =  I 

o 


-t  CO  ' 

D 


{wJ1(w)  +  UtJo(co)}2+  {coY1(co)  +  UtYo(co)}2 


dm 


co 


(C12 ) 


I  f(“)  d 

o 


CO 


Direct  numerical  evaluation  of  I,  is  difficult,  as  lim  f(co)  = 

t  co->0 

To  circumvent  this  difficulty,  the  Integral  was  evaluated  in 
two  parts . 

_3 

for  w  < 10  ,  the  following  situation  arises: 


coJ1  (co)  ->  0 


UTJo( ">  *  UT 


coY-^(co)  -*■  —0.6366 


W“> 

-fcD“2 


2Un 


->  1.0 


In  (co ) 


Thus 


f  (co) 


2  2UT  2 

w(U_  +  ( - -  In  hi  -  0.6366)  ) 

T  7T 


' 


■ 
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I 


W1 

/  f(to)  do) 
o 


1 


to=0 


where 


(0 


If  w  <  10  5, 
Thus  : 


2 


TT 


In  co  ~ 


0.6366 


to  (for  UT  =  00 )  <  -  7.34 


TT  dU) 

2Ut2  (1  +  w2) 


(Cl  3 ) 


to2  +  1  -  tO2 

-5 

If  to  is  taken  as  10  ,  the  error  introduced  by  the  approxima¬ 

tion  is  believed  to  be  negligible.  With  the  approximation, 
equation  C13  may  be  modified  and  integrated. 


I 


IT 

1  “ 

•> 

i 

o 

1 — 1 

II 

3 

IT 

1 

2 

2UtZ 

to 

O 

II 

3 

2U  2 

rp 

2  i 

—  In 

L7r 

0.6366 
to  -  — - 

UT  J 

to=10 


(0  =  0 


and 


I 


1 


2  -5 

-  In  (10  D) 

-IT 


0.6366 

Um 


(Cl  4 ) 


(to) 


UTJ0  <“)  } 


'+{u)Y1(m) 


UmY  (to)  } 
To 


dto 

_  (C15 

to 


Following  the  substitution  of  polynomial  approximations  for  the 
various  Bessel  functions  (page  369  in  Ref.  1),  the  integral 


, 


0  “W 


n  >  t  ’  f'.  ■  ■  -  ti'-  ’ 
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I,  was  evaluated  numerically  employing  a  nine  point  Gauss 
z2 

Quadrature  integration  formula.  In  an  attempt  to  optimize 
the  numerical  integration,  the  interval  size  (Aw)  was  increased 
as  w  increased. 

The  values  for  f(t)  calculated,  were  found  to  be 
comparable  to  those  presented  graphically  by  Ramey(ll). 


APPENDIX  II 


THE  COMPUTER  PROGRAM 
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A  computer  program  was  written  to  solve  the  four  equa¬ 
tions  developed  in  Appendix  I  simultaneously.  Upon  specifica¬ 
tion  of  time  and  steam  temperature;  depth,  quality,  casing 
temperature  and  f(t)  are  determined  by  trial  and  error  to 
within  specified  errors.  The  temperature  decrement  is  then 
adjusted  linearly  to  yield  an  increment  in  depth  approximately 
equal  to  the  specified  interval.  An  iterative  procedure  based 
on  successive  linear  adjustments  of  the  temperature  decrement 
is  applied  to  obtain  values  for  the  bottom  of  the  well.  Any 
numerical  integration  occurring  in  the  solution  was  performed 
using  a  nine  point  Gauss  Quadrature  formula.  Only  the  negative 
sample  points  and  their  corresponding  weighting  factors  are 
read  in,  the  other  values  are  generated  by  the  computer  with 
a  mirroring  procedure.  This  mirroring  routine  requires  the 
total  number  of  points  to  be  odd.  To  guard  against  endless 
loops  resulting  from  non-convergence,  both  subroutines  cal¬ 
culating  the  overall  heat  transfer  coefficient  contain  error- 
exits  . 

A  separate  program  was  employed  to  fit  the  tabulated 
data  for  the  heat  of  vaporization  and  CsL  to  polynomials  in 
temperature.  A  seventh  order  polynomial  yielded  a  maximum 
deviation  of  0.1%  between  the  calculated  and  tabulated  heat  of 
vaporization.  As  the  effect  of  C  T  on  the  overall  solution  is 

S  1j 

minor,  a  fourth  order  polynomial  fit  resulted  in  sufficient 
accuracy.  The  computer  program  is  limited  by  the  following 


'V5 

' 

■ 

;r  x  m  ci.vT 

■ 

... 
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restrictions : 

1)  50°C  <  T  <  374°C 

s 

2)  Pr*Gr  <  10^ 


3)  0.05  .<  y 
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I I -A.  FLOW  DIAGRAMS 


i)  Main  Program 


A 


B 


U 


. 
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E  F 


- - - — 
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+ 
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BO 


F 
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ii)  Subroutine  Coefficient  (single  casing) 


-  82 


A 


Is  ERROR 

<  ERROR  1 

yes 

-no 

Set  IXT  =  IXT  +  1 


Is  IXT 

v'  r 

yes 

<  5 

•  no 

Is  IN  =  2 


‘ 
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1 


The  flow  diagram  for  the  double  casing  situation  is  similar 
and  thus  not  included* 
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iii )  Subroutine  FT I ME 


Calculate  I, 

b 


>  Specify  an  interval  size  (DELT) 


Evaluate  the  incremental  area  of  the  interval  (Q) 
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II-R  FORTRAN  SOURCE  LISTS 


j)  MAIN  PROGRAM 


SIBFTC  ASTRI  NODECK 

COMMON  PI *UK  *DX  *  TR  *T IMEF  *LY*V ( 20 ) »H( 20 ) • ERT1 *  I T *R4*R6 * TC*UOfST* 

1AV I  SCI *AVISC2*AKl *AK2  *ACPl « ACP2 * ECON * ET ♦ ERC 1 »ERC2 *  A *B *C *D. A1 *A2  * 
2A3*AAA*BB*CC*E*F*FF*T3 

1  FORMAT ( 1 H  6F 1 3  «  6 ) 

2  FORMAT ( 1H  6E1 2«3 ) 

3  FORMAT ( I H  313) 

5  FORMAT  (1H  ^F14«,BS 

100  FORMAT ( 1HL16H  WELLBORE  AFTER  ?F10®2*5H  DAYS! 

1 01  FORMAT ( 1HL11H  INPUT  DATA/I 

102  FORMAT ( 1H  33HTEMPERATURE  LESS  THAN  LOWER  LIMIT?  3E14*6 ) 

103  FORMAT ( 1H  F 1 2 • 2  *  2 F 1 3 . 4 * E 1 4* 4 * F 1 0 • 4 » 2 F 1 3* 4  *  3 F 1 1 . 4 • I  5 ) 

110  FORMAT ( 1HK122H  DEPTH  STEAM  TEMP  HEAT  LOSS  CUM  HEAT  LO 

1SS  QUALITY  CASING  TEMP  EARTH  TEMP  F ( T  S  QV.  COEFF  @  F. 

2LOSS/ ) 

READ  (5*1 IACP1 *ACP2*AK1 * AK2 * AV I  SC  1  * AV I SC2 

READ  (5*1 )QA*QB*QUAL*WST*STEMP1 ?PA 

READ (5*1)  ECON *  EDENS *ECP*ETEMP1*ETEMP2*HEVAP 

READ( 5*1  )EM1 » EM2 * EM3 . EM4 * DEL T *  DEPTH 

READ  (5*1  )  F  I  LMC ♦ DTF I , CONDI » COND2 « COND3 * DX 

READ  ( 5*1  )  R 1 *R2 *R3*R4*R5*R6 

READ  ( 5  *2 )ERM1 *ERM2*ERM3  »ERM4*ERC1 *ERTl 

READ  (5*2) T0TIME*ADIST  *TMULT  *TMULT2*VGI *HRL IQ 

READ  (5*3)LM*NC*IWM 

READ  (5*5)  ( V ( L ) *  L= 1  *  LM ) 

READ  (5*5)  ( H ( L )  * L  =  1  *  LM ) 

WRITE (6*101) 

WR  I  TE ( 6*1 ) ACPI *ACP2*AK1 * AK2 * AV I  SCI  * AV I SC2 

WR  ITE( 6*1 )QA*QB*QUAL,WST*STEMPl ,PA 

WR  I  TE ( 6*1 ) ECON *  EDENS *ECP*ETEMPl,ETEMP2*HEVAP 

WR  ITE( 6*1 ) EMI *EM2*EM3*EM4*DELT, DEPTH 

WRITE* 6*1  )  F I LMC  *  DTF I  *  CONDI  *  COND2  *  COND3  *  DX 

WR  ITE(6*1 )Rl *R2*R3*R4*R5*R6 

WR  I  TE ( 6*2  )  ERMl *ERM2*ERM3*ERM4*FRC1*ERT1 

WR  I  TE ( 6  *  2 )T0TIME*ADIST*TMULT ♦ T MULT 2 ♦ VG I  *  HRL I Q 

WR  I  TE ( 6  *  3 ) LM  * NC  *  I WM 

WRITE  (6*5)  ( V ( L ) *  L= 1  *  LM ) 

WRITE  (6*5) (H(L) *L=1*LM) 

P I =3* 1415927 
L Y  =  2*LM-1 
L  T  =LY+ 1 
LS  =  LM+ 1 
DO  9  M  =  LS  » L Y 
LTM=LT-M 
V( M)s-V(LTM) 

9  H ( M ) =H ( LTM ) 

V ( LT )  =  0 • 

H ( LT  )  =  0 • 

C  CONSTANTS  FOR  HEAT  TRANSFER  COEFFICIENT 
HTLOSS=400.0 
U0=1.0 


■ 

, 


MW!  a  )  OA 

.  '  IT’  •* 
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FORTRAN  SOURCE  LIST  ASTRI 

T  3  =  ( STEMP1+ETEMP1 ) /2. 

A1=12./FI LMC/R1+AL0G(R2/R1 ) /CONDI 
A2  =  ALOG ( R4/R3 ) /COND2 
A  =  R  4  *  ( A1  +  A2 1/12® 

B=  R2*0 , 1713/ ( 1 •/EM1  +  R2*( 1./EM2-1* )/R3) /R4*0.  IE-07 
C= 12./ (R4*ALOG(R3/R2 ) ) 

D= ( R3-R? ) * ( (R3-R2)*28.9*PA/10«73)**2*32.2*0.75E+04 
CT  IME= 1 44 « *ECON/ ? EDENS*ECP*R4**2 ) *24 * 
CUO=R4/ECON/12® 

I F ( NC «  EQ® 1 )GO  TO  10 
A3=ALOG(R6/R5) /C0ND3 
AAA=R6*(A1+A2+A3)/12. 

BB=R4*B/R6 
CC=R4*C/R6 
E=12e/R6/ ALOG ( R5 /R4  5 

F=  R4*0 • 17 1 3/ ( l./EM3+R4*( 1./EM4-1. ) / R  5 ) /R6*0.1E-07 
FF  =  (R5'"R4  (  R5-R4  }*28«9*PA/10®73  }  **2*32 ® 2*0. 75E+04 

CT IME=CTIME*(P4/R6)**2 

CU0=R6/EC0N/12«  _  _ 

C  CONSTANTS  FOR  FRICTION  LOSS 

10  AA=3. 346313 
AB*0. 0414113 
AC=0»751548F-08 
AD  =  0 • 0 1 379448 1 

AE=0. 656444E-1 0  _ : _ 

A  T  =  647  #27 
Y=WST*6./R1/PI 

CZl=(WST/ (62.43*25.*Pl*Rl**2S  1**2/  C  64.4*777.97  ) 
CZ3=-0.2E~07/{ 5.04*1 .61*6.243*1 • 4^*R 1 **3 ) *Y**2 
C  ST-STEAM  TEMPERATURE 
C  Z TP=  D I  ST ANCE  DOWN  THE  WELLBORE 
C  OUAAa AVER AGE  QUALITY 
C  UO-OVERALL  HEAT  TRANSFER  COEFFICIENT 
C  T I MEF  =  F ( T  ) 

HEATIN=HRLIQ+QUAL*HEVAP 

TIME=DELT/TMULT-DELT*TMULT2 

11  T I ME=T I ME*TMULT+DELT*TMULT2  _ 

Z  T  P~0 . 

DZ  TP  =  AD  I S  T 

ST=STEMP1 

CHLOSS-O* 

HXLOSS=0. 

I  X  T  =  1 
DTF  =  DTF  I 
I  T  =  1 
INTC=0 
QUA2=QUAL 

V2S=(VGI*62.43)**2 
ST F=STEMP 1-460. 

ENTH2=HEVAP 
I W= I WM 

WRITE(6*100)TJME 
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fortran  source  list  astri 

WR I TE( 6*1 10) 

TR=TIME*CTIME 

12  UK=UO*CUO 
ER  R  =  UO 
CALL  FT  I  ME 

13  ET=ETEMP1+ZTP#ETEMP2 
CALL  COFFF  ( IXT ) 

I F ( ABS ( ERR /UO-1 « ) ®GT  ®ERM]  f  GO  TO  12 
IF ( I T , L T „ 2 ) GO  TO  21 

17  HTLOSS=Pl*ECON*( ( T Cl “FT 1 ) + ( TC-ET ) J/TIMFF 
IF(ABS(HL0SS/HTL0SS”l©)«LToERM4^G0  TO  22 
ZTP=ZTP~DZTP 
GO  TO  16 

22  CHLOSS=CHLOSS+HTLOSS*DZTP 
HXLOSS=HTLOSS 
21  ET  1  =ET 
I T  =  I T+l 
TC 1 =TC 

VI S=V2S  _ 

IF  <  ( IW.LT. IWM) .AND. ( INTC.LE.l ) ) GO  TO  30 
ETF=ET-460. 

TC  F  =  T  C-46  0. 

FLOSS  =  CHLOSS/WST /HEAT  I N 

WRITE(6*103) ZTP  *  STF  *  HXLOSS  *  CHLOSS  *QUA2*TCF.ETF*TIMEF.  UO*  FLOSS .  I  XT 
IW  =  0 

30  QUA1=QUA2 
I W= I W+l 

I F ( QUA 2 .LE.0.05JGO  TO  25 
IF (  INTC.GE.2 )GO  TO  25 
DTF=DTF*ADIST/DZTP 

14  STsST-DTF 
DTK=5.*DTF/9. 

ENTH1=ENTH2 
ST  F=ST-460. 

CPL=2.8317766-0. 0202964 39*$ TF+0. 81715712 E~0 4* STF**2~ 

10. 141 5405 1E-06*STF** 3+0. 921 29424E~1 0*STF**4 
EN TH2=0.1 591 4927E+04-0.1 040203 7E+02*STF+0. 8 1499083E-0l*STF**2- 
10. 3670295 8E-03*STF**3+0 . 96928323E-06*STF**4-0. 1 51 26743E~08*STF**5 
2+0 . 12907291E-11*STF**6-0.46713135E-15*STF**7 
SF  XP  =  0 . 

Xl=647.27-ST*5./9. 

X2  =X 1+DTK 

IF ( X2.GT.590. )GO  TO  24 
DO  15  1=1. LT 

XVa(Xl*(l.-V( I ) ) +X2* ( 1 ♦ +V ( I) ) )/2. 

IF ( I .EQ.LT ) XV=X2 

UE  =  XV*(  AA  +  AB*XV+AC*XV**3  +  AE*XV**4  »  /  (  (  1®+AD*XV)#(AT~XV)  ) 

DUE= ( AA+UE* ( l.-AT*AD  )+2»*( AB+UE*AD ) *XV+4 . *AC*XV**3+5. *AE*XV**4 ) / 
1( ( AT-XV ) * ( 1 . + AD*X V ) ) 

P=218.167*EXP(-UE*2.303) 

T  A  U= 1 • / ( AT-XV) 

G1 T=82.546*TAU-0.16246E+06*TAU**2 


, 

1  '■  T  l  I  •  T 
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FORTRAN  SOURCE  LIST  ASTRI 

G2T=0. 21828-0. 12697 E+06*TAU**2 

G3T=0. 3635 E-0 3-0.676 8 E~0 7*  (1 000. *TAU ) **24 

80  =  1,89-2641 .62*TAU*1  0.**(  8  08  70®  *TAU**2J 

B 1 =80* ( l.+BO*TAU*P*Gl T  + ( BO*T AU*P ) **3*G2T- « BQ*TAU*P ) ** 1 2*G3T ) 

VG  =  4«  55504/ ( P*TAU)+B1 

1 5  SF XP  =  SFXP-H( I )*2®3Q3*P*DUH/VG 
V2S=VG**2 

16  HLOSS=HTLOSS 

18  QUA2=(QUA1*I ENTH1+QUA1*CZ1*V15 )+CPL*DTF+ ( 1  • /777.97-HLOSS/WST )*DZTP 
1 )/ (ENTH2+QUA2*CZ1#V2S? 

QU AA= ( QUA 1+QUA2 ) /2. 

SZP=OZTP-SFXP/ (CZ3*QUAA**2*(®A+QB*QUAA) £*DTK 
IFUBS(SZP) •  LT, ERM2 ) GO  TO  19 

DE  R  =  1 . 0-SFXP / CZ  3*  < 2 • *QA+3. *QB*QUAA ) / ( QUAA**3* ( QA+QB*QUAA ) **2 ) * 

1 (HLOSS/WST-1 «/777®97 )  / ( 2 .* < ENTH2+QUA2*CZ 1*V2S > }*DTK 

DZTP=DZTP-SZP/OER 
GO  TO  18 

19  Z  TP  =  ZTP  +  DZTP 

___  IF  (INTC.EQ.2  )G0  TO  20 _ _ 

IF (ZTP.LT.DEPTH)GO  TO  13 

20  DELZTP=ZTP-DEPTH 
I N  TC  =  2 

IF (ABS(DELZTP) *LE.ERM3)GO  TO  13 
ST=ST+DTF 

V  ENTH2=ENTH1  _ 

ZTP=ZTP-DZTP 

DTF= ( 1,-DELZTP/DZTP) *DTF 
GO  TO  14 

24  WR ITE(6*102)ST*ZTP*X2 

25  IF (TIME.LT.TOTIME)GO  TO  11 

CALL  EXIT _ 

END 
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ii)  SUBROUTINE  FOR  THE  HEAT  TRANSFER  COEFFICIENT  (1  CASING) 

$  I  BFTC  COEFF1 

SUBROUTINE  COFFF  ( I  XT ) 

COMMON  PI ,UK,DX.TR.TIMEF»LY.V!20)  »H  ( 20) *  ERT  1 . IT  .R4.R6.TC.U0.ST  » 

1AV I  SCI *AVISC2  *AK1 ♦ AK2 . ACPI . ACP2 » ECON . ET ♦ ERC 1 . ERC2 t A .B .C , D . A1 . A2 • 
2A3.AAA.BB.CC.E.F.FF.T3 

10  FORMAT! 1HK40H  LIMIT  OF  K  CORRELATION  EXCEEDED  IN  C~1  *  El 4.6.2F10. 5 
1/) 

11  FORMAT ( 1HK23H  NO  CONVERGENCE  IN  Ol  *5E14.5/) 

I  N  =  I  T 

G=  R4*T I MEF/ 1 2  o 

1  T  2  =  ST-D I FF 1 
ERROR=l .0E+20 
I  X  T  =  1 

I N  = I N+ 1 

2  ERRORl =ERROR 
TM=!T2+T3)/2« 

DEL=T2-T3 

TMF=TM-460® 

AV  ISC  =  AVISC1+AVISC2*TMF 

AK=AK1+AK2*TMF 

CP=ACP1+ACP2*TMF 

DEN3=AK*AVISC*TM**3 

PRGR=D*DEL*CP/DEN3 

PAK=0.11*AK*PRGR**0.29 

IF ( (PRGR.LE.2000. ) . AND. ( I  XT . EQ • 1 )  ) PAK*AK 

3  DEN=A*B*( T2+T7)*(T2**2+T3**2)+A*C*PAK+1. 

UO=( l.-l . /DEN) /A 

T4  =  T  3-D  I FF2 

FT  =  T  4-ST+ECON*  (ST-ET  ) /(UO*G+ECON) 

ER  RORs ABS ! FT ) 

IF ( ERROR. LT.ERC1 ) GO  TO  4 

DPRGR  =  D*( DEL*ACP2-2.*CP-DEL*CP*( 3. /TM+AK 2/AK+AV ISC2/AV I  SC ) ) / 

1 ( DEN 3*2.7 

DPAKsPAK*! AK2/AK/2.+0.29*DPRGR/PRGR  5 

IF ( (PRGR.LE.200C. ) .AND. ( IXT.EO.l ) ) DPAK* AK2 / 2 . 

DEFT=1 ,-G*ECON*( ST~ET )*( B*( T2**2+2.*T2*T 3+3. *T3**2 )+C*DPAK) / ( DEN* 

1 (UO*G+ECON) ) **2 

T3  =T3-FT /DEFT  _ 

IF (ERROR. LT. ERRORl )GO  TO  2 
IXT=IXT+1 

IF  (  IXT.LE.6JGO  TO  2 

WR  ITE( 6*1 1 ) ERROR .ERRORl ,PRGRtPAK,AK 

4  DELQ=ECON*( T4-ET ) /TIMEF 
D I FF1=DELQ*A 1 

D I FF2=DELQ*A2 
IF (  IN.EQ.2IGO  TO  1 
6  TC  =  T4 

IF ! PRGR.GE.0.1E  +  07 ) WR I TE ( 6  *  1 0 ) PRGR ♦ P AK . AK 

RETURN 

END 
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iii)  SUBROUTINE  FOR  THE  HEAT  TRANSFER  COEFFICIENT  (2  CASINGS) 

SIBFTC  COEFF2 

SUBROUTINE  COEFF  (  IXT  ) 

COMMON  Pi  *  UK  * DX  *  TR  * T I MEF  Y ♦ V ( 20 ) »H( 20) ♦FRT1*IT*R4*R6,TC*U0,ET* 

1 A V  I  SC  1 ♦AVISC2.AKl?AK2*ACPl$ACP2eECON*ET»ERCl*ERC2*A»B*C.D*Al*A2* 
2  A3  ♦AAA»BB»CC»E«F?Fp!>T3 

10  FORMAT ( 1HK40H  LIMIT  OF  K  CORRELATION  EXCEEDED  IN  C~2  *2F14.6* 
14F10.5/ ) 

11  FORMAT ( 1HK23H  NO  CONVERGENCE  IN  C~2  *9E14.5/) 

I  N  =  I  T 

G=R4*T IMEF/12. 

GG=R6*T IMEF/12. 

1  T  2  =  $T-D I FF 1 
ERROR= 1 .0E+20 
I  X  T  =  1 

I  N  = I N+ 1 

2  ERROR1 =ERROR 
TM= ( T2+T3 ) /2. 

DEL=T2-T3 

TMF=TM-460.0 

AV ISC  =  AVI SC1+AVISC2*TMF  _ _ _ 

AK=AKl+AK2*TMF 

CP=ACP1+ACP2*TMF 

DEN3=AK*AVISC*TM**3 

PRGR=D*DEL*CP/DEN3 

PAK=0.11*AK*PRGR**0.29 

IF  (  (PRGR.LE.2000. )  .AND.J  IXT._E_0._1  IJ PAIC-AK _ _ _ _ _ _ 

3  FF 1  *  B*(T2+T3)*(T2**2+T3**2 ) 

DE  N 1 =  A*(FF1+  C*PAK ) +1 • 

UUO=(l. -l./DENl )/A 

T  4  =  T3-D I FF2 

T  6  =UUO* ( S  T-T4 ) *G/ECON+ET 

T  5  =  T  6+D I F  F3 _ _ 

D  E  L  A  =  A  8  $( T4-T  5 ) 

TMA=(T4+T5)/2. 

TMF=TMA-460.0 

AVISCA=AVISC1+AVISC2*TMF 

AKA=AK1+AK2*TMF 

CPA=ACP1+ACP2*TMF  _ _ 

DEN4=AKA*AVISCA*TMA**3 
PRGR2=FF*DELA*CPA/DENa 
PAK2-0.1 1 *AKA*PRGR2**0.  29 

IF ( ( PRGR2.LE. 2000. ). AND. ( IXT.EQ.l ) )PAK2=AKA 
FI =BB*FF1 /B+CC*PAK 

F2=F*( T4+T5)*(T4**2+T5**2 )+E*PAK2  _ 

DEN2=AAA*F1*F2+F1+F2 
UO= ( l.-(Fl+F2) /DEN2J/AAA 
FT=T6-$T+ECON*(ST~ET ) / ( UO*GG+ECON ) 

ER  ROR  =  ABS ( FT ) 

IF (ERROR. LT.ERC1 )GO  TO  4 

DPRGR  =  D*(  DEL*ACP2*-2.*CP“DEL*CP*(  3./TM+AK2/AK+AVISC2/AVISO  )  / 

1 ( D  EN3*2 . ) 

OPAK=PAK* ( AK2/AK/2.+0.29*DPRGR/PRGR ) 

IF ( { PRGR.LE.2000. ) .AND. ( IXT.EQ.l ) )DPAK*AK2/2. 


»  .*  -r'S* «  5T 
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FORTRAN  SOURCE  LIST  COFFF2 

IF (  (PRGR2»LE.2000. )  *  AND. (  I  XT  ®  EQ* 1 » )GO  TO  5 

DPRGR2=FF*(DELA*ACP2-2.*CPA-DELA*CPA*( 3. /TMA+AK2/ AKA+AV I SC2/AVISCA 
1 ) ) / ( 2.*DEN4) 

DPAK2=PAK2*( AK2/AKA/2.+0®29*DPRGR2/PRGR2 ) 
DPRGRS=DPRGR2+FF*2®*CPA/DEN4 
DPAK3=PAK2*( AK2/AKA/2.+0.29*DPRGR3/PRGR2 ) 

GO  TO  6 

5  DP AK2  =  AK2 /2  ® 

DPAK3=AK2/2. 

6  FF P  =  B*  (  T2**2-*-2#*T2#T3  +  3®*T3##2  ) 

ffpi=ffp+c*dpak 

FP 1 =FFP*BB/B+CC#DPAK 

DT  6  =  G* ( FFPl*( ST-T4 ) / DFNl **2-UUO ) /EOGN 

FP2=F*( 3.*T4**2+2.*T4*T5+T5**2 )+F*DPAK3+(F*(T4**2+2.*T4*T5+ 

1  3,*T5**2)+E*DPAK2)*DT6 

DUOT=( (F1+F2)*(AAA*( FP2*F1 +  FP1 #F2 ) +  FP1 +  FP2 ) /DFN2*«2- ( FP] +  FP2 ) /DEN2 
1 ) / AAA 

DEFT  =  DT6~ECON*GG*(ST**ET  )  /  (UO*GG+ECON  )**2*DUOT 
T  3  =T  3~F  T /DEFT 

IF ( ERROR. LT.ERRORl )GO  TO  2 
I  X  T= I XT  +  1 

IF ( I  XT • LE • 6 ) GO  TO  2 

WR I TE  C 6*11 ) ERROR* ERROR 1  * PRGR * PAK , AK ♦ PRGR2#PAK2* AKA 
4  DE.LQ  =  ECON*(T6-ET  )/TlMEF 
DI FF1=DELQ*A1 
01 FF2=DELQ*A2 
D I FF3=DELQ*A  3 
IF ( IN. EQ. 2 ) GO  TO  1 
TC  =  T6 

IF ( (PRGR.GT.1.0E+06) *0R. ( PRGR2 . GT . 1 . 0E+06 ) ) WR I TE( 6  *  10 ) PRGR ,PRGR2 * 
1AK  *PAKt AKA*PAK2 
RETURN 
END 
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iv  )  SUBROUTINE  FOR  THE  TIME  FUNCTION 

SIBFTC  TIME 

SUBROUTINE  FT  I  ME 

COMMON  Pl.UK.DX.TRtT I MEF  »l  Y  .  V  ( 20).H(20).ERT1*IT.R4*R6.TC*UO.ST. 

1 A V  I  SC 1 ♦AVISC2*AKl*AK2*ACPl*ACP2»ECON*ET»ERCl*ERC?*A*B*C*D*Al*A2* 
2 A3  ♦AAA*BB*CC*E*F*FF*T3 

SF=-PI/ (2. *UK**2*( 2. /PI *ALOG( 0.1 F-04) -0.6366/UK) ) 

Z2  =  0. 1F-04 
DELT=DX/TR 

1  Z  1  =Z  2 
Q=  0  . 

IF (BELT  .GT.0.001 ) GO  TO  2 

DELT=10.*DELT 

GO  TO  3 

2  OE LT= 1 • 2*DELT 

3  Z  2  =Z  2+DEL  T 
CON=0. 

DO  6  L  =  1 *  LY 

X* (Zl*( l.-V(L) )+Z2*(l.+VtL) ) ) /2. 

IF (X.GT.3. )GO  TO  4 
U*  X/3. 

FI *1.0-2. 2499997*0**2 
F2=X*( 0.5-0. 56249985*U**2 ) 

F3=2./Pl*ALOG( X/2. ) *F 1 +0 . 3674669 1+0. 605 59366*U**2 
F4=2./Pl*X*ALOG(X/2. ) *F2~0 . 6366 1 98+0 . 22 1 209 1 *U**2 
IF (U.LT.0.1E-02) GO  TO  5 

F1=F1+1.26?6208*U**4-0.3163866*U**6+0.0444479*U**R 
F2=F2+ ( 0. 2 1093 573-0.039542 89*U**2+0. 004433 19*U**4 ) *X*U**4 
F3 =F3-0.743S03P4*U**4+0. 25300 1 17*U**6~0. 04261 21 4*U**8 
F 4  =  F4+  (  2.  168270°-!  .  ^  1  64827*U**2+0.31  23951  *U**4  )  *!J**4 
IF (X.LT.0.1)  60  TO  5 

F 1 =F 1— ( 0. 3944E-02-0.21 F-03*U**2 ) *U**10 
F2=F2“ ( 0. 31^6!E-03-0.1 1 09E-04*U**2 ) *X*U**1 0 
F3=F3+( 0.427916E-02-0.24846E-03*U**2 )*U**10 
F4=F4-( 0. 400976E~01 “0.278 73 E~02*U**2 )*U**10 
GO  TO  5 

4  U  =  3  •  /  X 

FB =0.7978 8 456-0. 0000 0077*U~0. 00 55 2740*U**2“0. 0000 95 12*U**3+ 

10. 00 1 37237*U**4“0. 000728 05*U**5+0. 000 14476*U**6 
THETA  =  X~0.78939816~0. 04 1 663  97*U~0 . 0000  39  5 4*U** 2  +0 .00262577*11**3 
1-O.O00S41 2  5*U**4-0 . 00  0293  3  3*U**5  +  0  *  000 1 3  5  5  8  *U**6 
FBB=FB/X**0. 5 
FI =FBB*COS( THETA ) 

F3=FBB*SIN(THFTA) 

FB= 0.7978 8456+0. 000001 56*U+0. 01659667*U**2+0.0001 71 05*U**3- 
10. C0249cl i*u**4+0. 001 1365 3*U**5~0. 0002003 3*U**6 

TH ET A* X-2. 356 1 9449+0. 124996 12*U+0. 0000 565 0*U**2-0. 0063 7879*U**3 
1+0. 000 743 48*U**4+0. 00 0798 2 4*U**5-0. 000291 66*U**6 
FBB=FB/X**0.5 
F2=FBB*C05( THETA ) 

F4=FBB*X*SIN(THETA ) 

5  OEN=X*( (UK*F1+X*F2)**2+(UK*F3+F4)**2) 

XP=~TR*X**2 

IF ( A B S ( X P ) • G T . 4 0 • )  GO  TO  6 
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fortran  source  list  time 

Q::Q  +  OFLT*H(  L  )  *FXP  (  XP  )  /  (  2«*DEN  ) 

6  CONTINUE 
SF  =  SF+Q 

IF ( (0/SF*GT.FRTl ) .OR. (X*LT«2* ) )  GO  TO  1 
T I MEF= ( PI / <  2.*UK) ) #*2/SF~l • /UK 
RETURN 
END 
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II-C. 

NOMENCLATURE  OF  COMPUTER  INPUT 

ACPI 

Specific  heat  of  the  air  in  the  annulus  at  0°F 

ACP2 

(Btu/lbm°F) 

Specific  heat  change  with  temperature  (Btu/lb^  F  ) 

AKl 

Thermal  conductivity  of  the  air  at  0°F  (Btu/hr , ft . °F) 

AK2 

Thermal  conductivity  change  with  temperature 

(Btu/hr. ft.°F2) 

AVI  SC.1 

Viscosity  of  the  air  at  0°F  (lb  /hr. ft.) 

AVISC2 

Viscosity  change  with  temperature  ( lb^/hr . ft . °F) 

QA 

Coefficients  for  friction  factor  vs  quality  equation 

QS 

fG  =  QA  +  QB  »  quality 

QUAL 

Quality  of  the  inlet  steam  (fraction) 

WST 

Mass  flowrate  of  the  steam  (lb  /hr.) 

nr 

STEMPl 

Temperature  of  the  inlet  steam  (°R)  * 

PA 

Pressure  of  the  air  in  the  annulus  (psia) 

ECON 

Earth  conductivity  (Btu/hr. ft. °F) 

EDENS 

Earth  density  (lb^/ft2) 

ECP 

Earth  specific  heat  (Btu/lb  °F) 

m 

ETEMPl 

Earth  surface  temperature  (°R) 

E TEMP 2 

Geothermal  gradient  (°R/ft) 

IIEVAP 

Heat  of  vaporization  at  inlet  temperature  (Btu/lb^) 

EMI 

=  emissivity  of  outside  of  the  tubing 

EM2 

=  emissivity  of  inside  of  the  casing 

EM3 

£  ^  =  emissivity  of  outside  of  the  casing 

EM4 

=  emissivity  of  inside  of  the  second  casing 

i 
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. 
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DELT 

Basic  time  increment  (days) 

DEPTH 

Total  depth  of  the  well  (ft) 

FILMC 

Film  coefficient  for  the  inside  of  the  tubing 

(Btu/hr . ft. 2oF) 

DTFI 

Initial  steam  temperature  decrement  (°F) 

CONDI 

Thermal  conductivity  of  the  tubing  (Btu/hr . ft . °F) 

COND2 

Thermal  conductivity  of  the  casing  (Btu/hr . ft . °F) 

COND3 

Thermal  conductivity  of  the  second  casing 

(Btu/hr. ft . °F ) 

DX 

Initial  interval  for  evaluating  I 

t2 

Rl 

Inside  radius  of  the  tubing  (in) 

R2 

Outside  radius  of  the  tubing  (in) 

R3 

Inside  radius  of  the  casing  (in) 

R4 

Outside  radius  of  the  casing  (in) 

R5 

Inside  radius  of  the  second  casing  (in) 

R6 

Outside  radius  of  the  second  casing  (in) 

ERM1 

Specified  error  1  (see  main  flow  diagram) 

ERM2 

Specified  error  2  (see  main  flow  diagram) 

ERM3 

Specified  error  3  (see  main  flow  diagram) 

ERM4 

Specified  error  4  (see  main  flow  diagram) 

ERC1 

Specified  error  5  (see  subroutine  coefficient 

flow  diagram) 

ERTl 

Specified  error  6  (see  subroutine  FTIME  flow  diagram) 

TOTIME 

Maximum  injection  time  to  be  considered  (days) 

ADIST 

Approximate  depth  interval  to  be  considered  (ft) 

t 


TMULT 

coefficient  to  obtain  t  =  t* TMULT  +  DELT*TMULT2 

TMULT2 

coefficient  to  obtain  t  =  t* TMULT  +  DELT*TMULT2 

VGI 

Specific  volume  of  water  vapor  at  inlet  temperature 

(ft3/lb  ) 

HRLIQ 

Enthalpy  of  liquid  at  inlet  temperature  -  enthalpy 

LM 

of  liquid  at  earth  temperature  {Btu/lb^} 

Number  of  sample  points  for  Gauss  Quadrature  (neg.) 

(must  be  odd) 

NC 

Number  of  casings  (1  or  2) 

IWM 

Number  of  points  calculated  before  one  is 

printed  out 

V. 

l 

Sample  points  for  Gauss  Quadrature  (neg.  only) 

H  . 

l 

Weighting  factors  for  Gauss  Quadrature  for  above 

sample  points 

. 
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401  7  OJ  vD  3J<JMA2 


98 


The  results  of  one  computer  run  were  obtained  as  B.C.D. 
cards,  spaced  and  printed  on  the  line  printer,  and  subsequently 
included  in  this  section.  As  B.C.D.  card  output  is  restricted 
to  80  print  positions,  compared  to  the  line  printer's  132; 
f(t),  average  rate  of  heat  loss,  earth  temperature,  UQ  and  IXT, 
which  were  included  in  the  written  output,  were  not  punched 
out.  All  the  input  data  is  presented  as  a  print  out  of  the 
actual  computer  input,  thus  sections  II-B  and  II-C  are  required 
to  correlate  the  number  to  the  variable  and  its  units.  This 
section  also  contains  some  plots  of  friction  factor  vs  quality 
and  the  resulting  relationships.  An  attempt  to  fit  the  points 
of  Figure  15  to  a  polynomial  yielded  little  success,  attribut¬ 
able  to  the  odd  shape  of  the  curve.  As  a  result  straight  line 
approximations  were  made  for  all  cases,  rather  than  dividing 
the  curve  into  ranges. 


. noiJ-o  >4  a  li  r.l  F'  )LiJ.  oni 


/ 


;  ‘  h  ,  Ioj  IsinonYlcq  t  o  5f  vac/pj  :o 
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FIG.  14 .  TWO  PHASE  FRICTION  FACTOR  VERSUS 
VAPOR  PHASE  REYNOLDS  NUMBER 
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III-A.  TWO  PHASE  FRICTION  FACTOR 


Figure  15  Data: 

W  .  =  steam  flowrate  =  3000  lb  /hr 

st  m 


Injection 

pressure 

=  500  psia 

Tubing  diameter  (inside)  - 

2.441  in 

From  Keenan  and  Keyes (8)  obtain  the  values  for 

saturated 

viscosity 

and  liquid  density. 

^  G 

=  5.45  x 

10  *™2  lb  /hr- ft 
m 

PL 

=  50.7  lb 

/ft3 

nr 

Thus  : 

Y 

2Wst 

2  x  3000 

a 

’Yg 

,  2.441.  _  .  c 
n  2  x  12  5,45 

x  10~2 

=  3.45  x 

105 

v 

W  . 
st 

3000 

=  92.5 

Yb 

2"rlPL 

11  < 2 12 4 1  >  50  *  7 

and 

=  3. 45  x 

105  y 

Y 

=  92.5  (1 

-  y) 

With  the  use  of  Figure  14,  values  for  the  two  phase  friction 
factor  may  be  determined  for  various  qualities 


. 


FIG.  15.  TWO  PHASE  FRICTION  FACTOR  VERSUS  QUALITY 
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Table  I 


y 

’’eg  x  10~4 

Y 

f  G 

0.05 

1.73 

87.8 

0.043 

0.10 

3.45 

83.3 

0.044 

0.15 

5.16 

78.6 

0.0  31 

0.20 

6 . 90 

74.0 

0.026 

0.30 

10.4 

64.8 

0.023 

0.40 

13.8 

55.5 

0.018 

0.50 

17.3 

46.3 

0.016 

0.60 

20.7 

37.0 

0.014 

0.70 

24.2 

27.8 

0.012 

0.80 

27.6 

18.5 

0.010 

0.90 

31.1 

9.25 

0.008 

1.00 

34.5 

0 

0.0035 

Figure  17  Data: 

W  .  =  10 , 500  lb  /hr 

st  m 

Average  pressure  -  700  psia 

\i  G  =  0.0  8  lb^/hr.ft.  (from  Perry  (10)) 

pT  =  48.8  lb  /ft3 

L  m 


Part  i) 


Tubing  diameter 


Y 

a 


=  1.995  in 


2W 


s  t 


"rlUG 


2  x  10,500 


1  995 

7T  (-irJLld.)  o.0  8 


1.005  x  106 


2  x  12 


oe.o 
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2TrrlPL 


10,500 


7T  ( 


1.995 

12 


)  48.8 


411 


Thus  for  Part  i)  have: 


Req  =  1.005  x  106  y 

Y  =  411  (1  -  y) 


Part  ii) 

Tubing  diameter  =  2.441  in 

Reg  =  8.21  x  105  y 

Y  =  337  (1  -  y) 


Part  iii) 

Tubing  diameter 

Y 


2.992  in 
6.69  x  105  y 
274  (1  -  y) 


The  y  values  involved  are  out  of  the  range  of  Figure 
14,  thus  Figure  16  was  obtained  and.  the  curves  extrapolated. 
The  values  so  obtained  were  plotted  on  Figure  14  and  a  curve 
was  then  sketched  in. 
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X 

Qa 

rtf 

U 

tn 

CD 

I 

w 

OS 

jC 

-P 

c 

o 

t>1 

-P 

•H 

i — I 

rtf 

P 

D1 

•P 

CO 

a 

-H 

<T3 

tn 

rtf 

n3 

05 

-P 

■P 

O 

l — I 

CU 

0) 

u 

rtf 

CO 

-P 

U 

rtf 

Cu 

Q) 

0) 

U 

X 

-P 


rtf 

B 

O 

U 

4-1 

CO 

u 

O 

-P 

O 

rtf 

4-i 

C 

O 

•iH 

-P 

o 

•H 

U 

4-1 

Q) 

XI 

En 


U 

CD 

-P 

4-' 

rtf 

O 

CO 

<D 

X 

Eh 


of  points  was  such  that  one  straight  line  could  reasonably  fit  all  three  parts. 
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I  I  I  —  B  COMPUTER  INPUT  DATA 


FIG  2  AND  3 

BASIC  INPUT  DATA  (ALSO  INPUT  RESULTING  IN  SAMPLE  OUTPUT) 


0.239 
0.0285 
1 . 0 
0.8 

2000.0 
1.2205 
0.500E-01 
0.  100E  +  04 
5  1  4 

-0.96816024 
0.08127439 


0.000018 

-0.0218 

110.0 

0.8 

0.5 

1.4375 
0. 100E-01 
0.250E+02 


0.013 
1  .00 
0.197 
0.8 
25.0 
2.9945 
o.iooE+no 
0  « 1 00E  +  02 


0*000026 
5000.0 
535.0 
0  ®  8 
25.0 
3.3125 
0. in^E-03 
0.0 


0.04 
927.0 
0.0110 
3.65 
25.0 
5.025 
0  ® 1 onE~n2 
9.278E-01 


0.000056 


14.7 
755.0 
6000*0 
0.00001 
5.375 
0. 100E-04 
4®064E+02 


-0.83603111 

0.18064816 


-0*61337143 

0*26061070 


*0*32425342 

0.31234708 


-0.0 

0.33023936 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


CHANGES;  QB  TO  -0.0225  AND  WST  TO  3000.0  ON  CARD  2 


FIG  4 

BASIC  INPUT  DATA 


0.239 
0.0285 
1.4 
0.8 
2  J00 • 0 

1.2205 
0.300E-O1 
0. 500E+01 
5  1  4 

-0.96816024 
0.08127439 


0.000018 
.  -0.0218 
140.0 
0.8 
0.3 

1.4375 
0. 100E-01 
0.250E+02 


0.013 
0.90 
0.23 
0.8 
25.0 
1.979 
0.100E+OO 
0. 100E  +  01 


0.000026 
4850.0 
518.0 
0.8 
25.0 
2.250 
0. 100E-03 
0.  100E  +  01 


0.040 
862.0 
0.0263 
2.0 
25.0 
5.025 
o.looE-02 
1 .823E+00 


0*000056 


14.7 

824.1 

1637.0 

0.00001 
5.375 
0. 100E-04 
3.500E+03 


-0.83603111 

0.18064816 


-0.61337143 

0.26061070 


-0.32425342 

0.31234708 


-0.0 

0.33023936 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


CHANGES?  QUAL  ON  CARD  2  TO  0.8 


FIG  5.6  AND  7 


BASIC  INPUT  DATA 


0.240 

0.0285 

0.4 

0.8 

2000.0 
0.9975 
0. 500E-01 
0. 100E+05 
5  1  4 

-0.96816024 
0.08127439 


0.000023 

-0.0213 

140.0 

0.8 

0.5 

1.1875 

0.100E-01 

0.250E+02 


0.0133 
0.80 
0.23 
0.8 
25.0 
2.522 
0. 100E  +  00 
0. 100E+02 


0.000021 
10500.0 
520.0 
0.8 
25.0 
2.750 
0. 100E-03 
0.0 


0.0395 
1020.0 
0.0108 
1.00 
25.0 
5.025 
0.100E-02 
3.868E-01 


0.000050 


14.7  2 

624.2  3 

3250.0  4 

0.00001  5 

5.375  6 

O.100E-04  7 

5.342E+02  8 

9 


-0 .83603111  -0.61337143  -0.32425342  ~0.0  10 

0.18064816  0.26061070  0.31234708  0.33023936  11 


' 


■ 


:  :  •  i  . 


e?*TSf8Q. 
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CHANGES;  R 1 

TO 

1 

.1250 

AND 

R2 

TO 

1 

.3125 

ON 

CARD 

6 

R  1 

TO 

1 

.2205 

AND 

R2 

TO 

1 

.4375 

ON 

CARD 

6 

R 1 

TO 

1 

.  3750 

AND 

R  2 

TO 

1 

.6250 

ON 

CARD 

6 

R1 

TO 

1 

.4960 

and 

R2 

TO 

1 

,7500 

ON 

CARD 

6 

FIG  8  AND 

9 

BASIC  INPUT  DATA 

0.240 

0.000023 

0.0133 

0  ®  00002 1  0.0395 

0.000050 

0.0285 

-0.0218 

0.80 

8000.0  1047.0 

14.7 

2 

0.4 

140.0 

0.23 

520.0  0.0108 

574,7 

3 

0.8 

2000.0 

0.8 

0.300 

0.8 

25.0 

0.8  1.00 

25.0  25,0 

3250,0 

o.ooool 

4 

5 

1.22055 

1.4375 

2.045 

2.250  5.025 

5,375 

6 

0.500E-'M 

0.10OE-01  0 

.  1O0E+0O 

o.looE~02  <">  •  1  OOE-n  1  n. 

1  ooE-04 

7 

0. 100E+04 

0. 100E+03  0 

•  1 00E  +  02 

0.0  3.012E-01  5. 

707E+02 

8 

5  1  1 

9 

-0.96816024 

-0.83603111 

-0.61337143  -0.32425342  -0.0 

10 

0.08127439 

0.18064816 

0.26061070  0.31234708  0.33023936 

11 

CHANGES; 

WST  ON  CARD  2 

TO  10500* 

12000,  14000*  18000 

FIG  10*11,12  AND  13 


BASIC  INPUT  DATA 


0.240 
0. 0285 
0.4 
0.8 

2000.0 
0.9975 
0. 500E-01 
). 100E+05 
5  14 

-0.96816024 
0.08127439 


0.000023 

-0.0218 

140.0 

0.8 

0.5 

1.1875 
0. 100E-01 
0.25CE+02 


0.0133 
0.80 
0.23 
0.8 
25.0 
2.045 
0. 100E  +  00 
0. 100E+02 


0.000021 
10500.0 
520.0 
0.8 
25.0 
2.250 
0. 100E-03 
0.0 


0.0395 
1020.0 
0.0108 
1.00 
25.0 
5.025 
0.100E-O2 
3.868E-01 


0 .000050 


14.7 

624.2 

3250.0 

0.00001 
5.375 
0. 100E-04 
5.342E+02 


-0.83603111 

0.18064816 


-0.61337143 

0.26061070 


-0.32425342 

0.31234708 


-0.0 

0.33023936 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


CHANGES: 


STEMP1  ON  CARD  2  TO  1040.0 

HEVAP  ON  CARD  3  TO  588.4 

VG I  ON  CARD  8  TO  3.217E-01 

HRLIQ  ON  CARD  8  TO  5.609E+02 


*  1060.0 
*  548.5 

,  2.668E-01 
,  5.890E+02 


,  1080.0  *  1100. 0 
»  503.6  ,  452.0 

,  2.201E-01  ,  1.798E-01 
*  6.187E+02  ,  6.506E+02 


NOTE;  ONE  COLUMN  REPRESENTS  ONE  DIFFERENT  RUN 


THE  CHANGES  FOR  FIG  13  ARE  STATED  ON  THE  GRAPH 


•  '  ••  • 
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Ill-c  SAMPLE  COMPUTER 

OUTPUT 

(  for 

FIG  2  AND 

3  ) 

WELLBORE 

AFTER 

3,65  DAYS 

DEPTH 

STEAM  TEMP 

CUM  HEAT 

LOSS 

QUALITY 

CASING  TEMP 

F.  LOSS 

(  FT  ) 

(  F) 

( BTU/HR ) 

(  F) 

0.00 

467  «  '"'000 

OeOOOOE- 

•38 

1 *0000 

279.1233 

0 .0000 

1 67.06 

466.0918 

0.8600E 

05 

0.9774 

279.3994 

0.0148 

266.52 

465.5360 

0.1369E 

06 

0.9641 

279.5535 

0.0236 

366.00 

464.9684 

0.1876E 

06 

0*9509 

279.6992 

0*0323 

465.54 

464.3894 

0.2381E 

06 

0  ®  9  078 

279.8375 

0.0410 

565.12 

463.7999 

0.2888E 

06 

0.9247 

279,9689 

0,0496 

664.73 

463.2004 

0®  3383E 

06 

0.91  18 

280.0936 

0.0583 

764.37 

462.5920 

0. 388  IE 

06 

0,8989 

280.2124 

0.0668 

864.04 

461.9749 

0.4377E 

06 

0.8861 

280.3257 

0*0754 

963.74 

461.3501 

0.4870E 

06 

0.8735 

280,4342 

0.0839 

1063.47 

460.7178 

0.5360E 

06 

0.8609 

280.5383 

0.0923 

1163.23 

460.0790 

0.5849E 

06 

0.8483 

280.6384 

0*1007 

1263.02 

459.4342 

0.6335E 

06 

0.8059 

280.7350 

0.1091 

1362.82 

458.7841 

0.6819E 

06 

0.8236 

280.8287 

0.1174 

1462.66 

458.1292 

0.7300E 

06 

0.8113 

280.9198 

0.1257 

1562.51 

457.4701 

0.7779E 

06 

0.7092 

281.0090 

0.1340 

1662.39 

456.8075 

0.8255E 

06 

0.7871 

281.0967 

0. 1422 

1762.29 

456.1419 

0.8729E 

06 

0.7752 

281,1831 

0,1503 

1862.22 

455.4738 

0.9200E 

06 

0.7633 

281.2687 

0.1584 

1962.16 

454.8039 

0.9669E 

06 

0.7515 

281.3543 

0.1665 

2062.12 

454.1327 

0.1014E 

07 

0.7398 

2  e 1 .4396 

0.1745 

2162.11 

453.4608 

0. 1060E 

07 

0.7282 

281.5258 

0.1825 

2262.11 

452.7886 

0. 1 106E 

07 

0.7166 

281.6128 

0.1905 

2362.13 

452.1169 

0.1152E 

07 

0.7052 

281.7012 

0.1984 

2462.18 

451.4458 

0.1198E 

07 

0 • 6°38 

281.7913 

0.2062 

2562.24 

450.7762 

0.1243E 

07 

0.6826 

281.8834 

0.2141 

2662.33 

450.1084 

0. 1288E 

07 

0.6714 

281.9779 

0.2218 

2762.43 

449.4431 

0.1333E 

07 

0.6603 

282.0753 

0.2296 

2862.55 

448.7806 

0.1378E 

07 

0.6493 

282.1758 

0.2373 

2962.69 

448.1215 

0.1422E 

07 

0.6383 

282.2799 

0.2449 

3062.85 

447.4663 

0.1466E 

07 

0.6275 

282.3877 

0.2525 

3163.02 

446.8155 

0.1510E 

07 

0.6167 

282.4998 

0.2601 

3263.22 

446.1693 

0.1554E 

07 

0.6060 

282.6164 

0.2676 

3363.44 

445.5285 

0.1598E 

07 

0.5954 

282.7377 

0.2751 

3463.68 

444.8935 

0.1641E 

07 

0.5849 

282.8643 

0.2825 

3563.93 

444.2646 

0.1684E 

07 

0.5744 

282.9961 

0.2899 

3664.21 

443.6424 

0.1726E 

07 

0.5641 

283.1336 

0.2973 

3764.51 

443.0272 

0  •  1 76QE 

07 

0.5^38 

283.2773 

0.3046 

3864.82 

442.4195 

0  •  1 8  1 1 E 

07 

0.5435 

283.4270 

0.3119 

3965.17 

441.8196 

0.1853E 

07 

0.5334 

283.5834 

0.3191 

4065.53 

441.2280 

0.1895E 

07 

0.5233 

283.7466 

0.3263 

4165.91 

440.6450 

0. 1937E 

07 

0.5133 

283.9168 

0.3335 

4266.32 

440.0712 

0. 1978E 

07 

0.5033 

284.0942 

0*3406 

4366.74 

439.5067 

0.2019E 

07 

0.4934 

284.2791 

0.3477 

4467.20 

438.9519 

0.2060E 

07 

0.4836 

284.4717 

0. 3547 

4567.66 

438.4074 

0.2101E 

07 

0.4738 

284.6722 

0.3618 

4668.17 

437.8731 

0.2141E 

07 

0.4642 

284.8809 

0.3687 

, 

*  - 

3PST8.0 

30C  i^.O 

w.o 

vsrr.^e# 

:  •  ) 

0  1  «  • 

s  /  .  )?fi 

.  %  Is  ?  * 

A>8r.%44 

'IS  £.844 

'  i  ,  0 

34<3'i  *>  4 

-  .0 
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WELLBORF  AFTER  3.65  DAYS 


DEPTH 

STEAM  TEMP 

CUM  HEAT 

LOSS 

Quality 

CASING  TEMP 

F.  LOSS 

(  FT) 

(  F) 

( BTU/HR  ) 

(  F) 

4768.70 

437.3496 

0  •  2 1  8  1 E 

07 

0.4545 

285.0978 

0.3757 

4860.24 

436.8372 

0  @  22  2 1 E 

07 

0.4449 

285.3232 

0.3825 

4969.82 

436.3360 

0.2261E 

07 

0.4354 

285.5571 

0.3894 

5070.42 

435.8465 

0.2301E 

07 

0.4260 

285.7048 

0.3962 

5171.06 

435.3688 

0.2340E 

07 

0.4166 

285.9572 

0.4030 

5271.71 

434.9032 

0.2380E 

07 

0.4072 

286.2184 

0.4098 

5372.40 

434,4499 

0  ®  24 1 9E 

07 

0.3979 

286.4888 

0.4165 

5473.12 

434.0090 

0  ®  2457E 

07 

0.3886 

286.7686 

0.4232 

5573.87 

433.580° 

0.2496E 

07 

0.3794 

287.0578 

0.4298 

5674.65 

433.1656 

0.2S34E 

07 

0.3703 

287.3562 

0.4364 

5775.48 

432.7633 

0.2573E 

07 

0.3612 

287.66*4 

0.4430 

5876.32 

432.3741 

0.2611E 

07 

0.3521 

287.9818 

0.4496 

5977.22 

431.9981 

0.2648E 

07 

0.3431 

288.3092 

0.4561 

5999.99 

431.9151 

0.2657F 

07 

0.3411 

288.3843 

0 • *5  7  5 

WELLBORE 

AFTER 

36.50  DAYS 

0.00 

467.0000 

0.0000E- 

■38 

1.0000 

318.0724 

0 • 0000 

167.25 

466.0928 

0.7083E 

05 

0.9814 

318.0571 

0.0122 

266.74 

465.5395 

0.U27E 

06 

0.9705 

318.0373 

0.0194 

366.26 

464,9753 

0.1544E 

06 

0.9597 

318.00°0 

0.0266 

465.81 

464.4002 

0. 1959E 

06 

0.9489 

317.9723 

0.0337 

565.38 

463.8149 

0.2372E 

06 

0.9382 

317.9276 

0.0409 

664.98 

463.2197 

0.27S4E 

06 

0.9276 

317.8757 

0.0479 

764.60 

462.6151 

0.3193E 

06 

0.9170 

317.8164 

0.0550 

864.24 

462.0014 

0. 3601 F 

06 

0.9065 

317.7506 

0.0620 

963.90 

461.3790 

0.4006E 

06 

0.8961 

317.6782 

0.0690 

1063.59 

460.7483 

0 • 44 1 OF 

06 

0.8858 

317.5999 

0.0759 

1163.29 

460.1097 

0.4812E 

06 

0.8755 

317.5159 

0.0829 

1263.00 

459.4635 

0. 52 1  IE 

06 

0.8653 

317.4265 

0.0897 

1362.73 

458.8101 

0.5609E 

06 

0.8552 

317.3320 

0.0966 

1462.48 

458.1499 

0.6004E 

06 

0.8452 

317.2330 

0.1034 

1562.25 

457.4831 

0.6398E 

06 

0.8352 

317.1295 

0.1102 

1662.03 

456.8102 

0.6789E 

06 

0.8253 

317,0220 

0.1169 

1761.82 

456.1314 

0.7179E 

06 

0.8155 

316.9107 

0.1236 

1861.63 

455.4471 

0.7566E 

06 

0.805R 

316.7958 

0.1303 

1961.44 

454.7576 

0.7952E 

06 

0.7961 

316.6778 

0.136° 

2061.28 

454.0632 

0.8335E 

06 

0.7866 

316.5568 

0.1435 

2161.12 

453.3643 

0.8716E 

06 

0.7771 

316.4332 

0.1501 

2260.98 

452.6611 

0.9095E 

06 

0.7676 

316.3073 

0.1 566 

2360.84 

451.9539 

0.9472E 

06 

0.7583 

316.1793 

0.  1631 

2460.72 

451.2432 

0.9847E 

06 

0.7490 

316.0495 

0.1696 

2560.61 

450.5291 

0. 1022E 

07 

0.7398 

315.9182 

0.1760 

2660.51 

449.8119 

0. 1059E 

07 

0.7307 

315.7855 

0.1824 

2760.42 

449.0920 

0. 1096E 

07 

0.7216 

315.6516 

0.1987 

2860.34 

448.3696 

0.1 1 33E 

07 

0.7127 

315.5173 

0.1°50 

2960.27 

447.6450 

0.11 69E 

07 

0.7038 

315.3821 

0.2°17 

3060.21 

446.9185 

0. 1 205E 

07 

0 • 6°50 

715.2466 

0.207* 
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WELLBORE  AFTER 


36.50  DAYS 


DEPTH 

STEAM  TEMP 

CUM  HEAT 

LOSS 

QUALITY 

CASING  TEMP 

F.  LOSS 

(  FT  ) 

(  F) 

«  BTU/HR ) 

(  F) 

3160.16 

446.1903 

0.1241F 

07 

0,6862 

315.1113 

0.2137 

3260.12 

445.4608 

0.1277E 

07 

0*6775 

314.9759 

0.2199 

3360.09 

444.7303 

0.1313E 

07 

0.6689 

314.8409 

0.2260 

3460.08 

443.9989 

0.1348E 

07 

0.6604 

314.7065 

0.2321 

3560.05 

443.2671 

0. 1383E 

07 

0.6520 

314,5731 

0.2382 

3660.05 

442.5348 

0.141 8E 

07 

0.6436 

314.4406 

0,2442 

3760.06 

441.8026 

0.1453E 

07 

0.6353 

314.2577 

0,2502 

3860.07 

441.0707 

0. 1487E 

07 

0.6270 

314. 12e3 

0.2561 

3960.09 

440.3393 

0.1522E 

07 

0.6189 

314.0005 

0.2620 

4060.13 

439.6087 

0.1556E 

07 

0.6108 

313,8746 

0.2679 

4160.17 

438.8792 

0.1590E 

07 

0.6027 

713.7509 

0.2737 

4260.22 

438.1509 

0.1623E 

07 

0.594R 

313,6294 

0.2795 

4360.28 

437.4242 

0.1657E 

07 

0.5869 

313,5105 

0.2853 

4460.34 

436.6993 

0.1690F 

07 

0.5791 

313.3943 

0.2910 

4560.42 

435.9764 

0. 1723E 

07 

0.5713 

313.2809 

0.2967 

4660.50 

435.2559 

0.1756E 

07 

0.5636 

313.1707 

0.3024 

4760.60 

434.5378 

0. 1789E 

07 

0.5560 

313.0636 

0,3080 

4860.70 

433.8226 

0. 1821E 

07 

0.5485 

312.9601 

0.3136 

4960.81 

433.1105 

0.1853E 

07 

0.5410 

312.8602 

0.3191 

5060.93 

432.4016 

0.1885E 

07 

0.5336 

312.7642 

0.3247 

5  161. .06 

431.6962 

0.1917E 

07 

0.5262 

312.6721 

0.3301 

5261.20 

430.9946 

0.1949F 

07 

0.5189 

312.5841 

0.7356 

5361.35 

430.2969 

0.19R0E 

07 

0.5117 

312.5005 

0.7410 

5461.51 

429.6036 

0.201  IE 

07 

0.5045 

312.4215 

0.3464 

5561.68 

428.9146 

0.2042E 

07 

0.4974 

312.3471 

0.3517 

5661.85 

428.2305 

0.2073E 

07 

0.4904 

312.2776 

0.3570 

5762.03 

427.5511 

0.2104E 

07 

0.4834 

312.2130 

0.3623 

5862.23 

426.8769 

0.2134E 

07 

0.4764 

312.1536 

0.3675 

5962.44 

426.2081 

0.2164E 

07 

0.4696 

312.0995 

0.3727 

6000.01 

425.9589 

0.2175E 

07 

0.4670 

312.0807 

0.3746 

WELLBORE  AFTER  365.00  DAYS 


0.00 

467.0000 

o.ooooe- 

-38 

1.0000 

344.7480 

o.oo^o 

167.38 

466.0936 

0.5985E 

05 

0.9843 

344.1456 

0.0103 

266.91 

465.5422 

0.9521E 

05 

0.9751 

344,0148 

0.0164 

366.46 

464.9805 

0. 1304E 

06 

0.9660 

343.8755 

0.0225 

466.02 

464.4086 

0.1655E 

06 

0.9569 

343.7279 

0.0285 

565.60 

463.8268 

0.2004E 

06 

0.9479 

343.5721 

0.0345 

665.20 

463.2354 

0.2351E 

06 

0.9390 

343.4085 

0.0405 

764.82 

462.6345 

0.2697E 

06 

0.9301 

343.2371 

0.0464 

864.45 

462.0245 

0.3041E 

06 

0.9212 

343.0583 

0.0524 

964.09 

461.4055 

0.3383E 

06 

0.9125 

342.8724 

0.0583 

1063.75 

460.7777 

0.3724E 

06 

0.9038 

342.6797 

0.0641 

1163.42 

460.1414 

0.4063E 

06 

0.8951 

342.4802 

0.0700 

1263.10 

459.4968 

0.4401E 

06 

0.8866 

342.2741 

0.07SP 

1362.80 

458. 8440 

0.4736E 

06 

0.8780 

342.0618 

0.0816 

1462.50 

458.1833 

0. 5070E 

06 

0.8696 

341  .84^3 

0.0873 

Y  A  0  0  *  • 

A  yt  ■ 

. 

P  P  '=.  A 

Vi 

- 

J4C- 

'  F  'A* 

a*.  ,0 

;  >S*;  !*fi. 

•’  V  >r  . 

*TP*F£A 

<V  •  S  F  A 

■V  .  'P? 

:•  ♦’.  >* 

t  .  ♦' 

WELLBORE  AFTER 


365.00  DAY A 


DEPTH  STEAM  TEMP  CUM  HEAT  LOSS  QUALITY  CASING  TEMP  F.  LOSS 


(  FT  ) 

(  F) 

( fITM/HR  ) 

1562.22 

457.5148 

0.5403E 

06 

1661.94 

456.8387 

0  ®  5 73 3  E 

06 

1761.67 

456.1552 

0.6062E 

06 

1861.42 

455.4645 

0.6389E 

06 

1961.17 

454.7667 

0.6714E 

06 

2060.93 

454.0620 

0.7038E 

06 

2160.69 

453.3505 

0.7360E 

06 

2260.47 

452.6322 

0.7680E 

06 

2360.25 

451.9076 

0.7998E 

06 

2460.04 

451.1766 

0.8314E 

06 

2559.83 

450.4395 

0.8629E 

06 

2659.63 

449.6962 

0.8942E 

06 

2759.43 

448.9470 

0.9253E 

06 

2859.24 

448.1920 

0.9562E 

06 

2959.06 

447.4313 

0.9869E 

06 

3058.88 

446.6650 

0.1 017E 

07 

3158.70 

445.8932 

0. 1048E 

07 

3258.53 

445.1160 

0. 1078E 

07 

3358.36 

444.3335 

0. 1 108E 

07 

3458.20 

443.5458 

0.1 138E 

07 

3558,04 

442.7530 

0.1 168E 

07 

3657.88 

441.9552 

0 • 1 1 97E 

07 

3757.73 

441 .1524 

0.1226E 

07 

3857.58 

440.3448 

0.1255E 

07 

3957.44 

439.5324 

0.1284E 

07 

4057.29 

438.7153 

0.1313E 

07 

41 57.15 

437.8935 

0.1341E 

07 

4257.02 

437.0671 

0.1370E 

07 

4356.88 

436.2362 

0.1398E 

07 

4456.75 

435.4008 

0.1426E 

07 

4556.61 

434.5610 

0.1454E 

07 

4656.48 

433.7167 

0. 1481E 

O7 

4756.36 

432.8681 

0.1508E 

07 

4856.23 

432.0152 

0.1536E 

07 

4956.10 

431.1580 

0. 1 563E 

O7 

5055.98 

430.2964 

0.1589E 

O7 

5155.86 

429.4307 

0. 1 61 6E 

07 

5255.73 

428.5606 

0.1642E 

O7 

5355.61 

427.6864 

0.1669E 

O7 

5455.49 

426.8079 

0.1695E 

07 

5555.37 

425.9251 

0.1720E 

07 

5655.25 

425.0381 

0. 1746E 

07 

5755.13 

424.1468 

0.1772E 

07 

5855.00 

423.2512 

0.1797E 

07 

5954,88 

422.3513 

0.1  8  2  2E 

07 

6000.02 

421.9435 

0. 1 833F 

O7 

<  F) 


0.8612 

341.6189 

0 

.0930 

0.8529 

341.3889 

0 

.0987 

0.8446 

341.1533 

0 

.1044 

0.8364 

340.9123 

0 

.1100 

0.8283 

340.6660 

0 

9  1 1  56 

0.8203 

340.4149 

0 

.1212 

0.8123 

340.1588 

0 

.1267 

0.8043 

339.8979 

0 

.1322 

0.7965 

339.6326 

0 

.1377 

0.7887 

339,3629 

0 

.1432 

0.7809 

339.0889 

0 

.  I486 

0.7733 

338.8108 

0 

.1540 

0.7657 

338.5286 

0 

,1593 

0.7582 

338.2424 

0 

.1647 

0.7507 

337.9528 

0 

.1700 

0.7433 

337.6594 

0 

.1752 

0.7360 

337.3626 

0 

.1804 

0.7287 

337.0622 

0 

.1356 

0.7215 

336.7309 

0 

.1908 

0.7144 

336.4243 

0 

.1959 

0.7073 

336.1147 

0 

.2011 

0.7003 

335.8021 

0 

.2061 

0.6934 

335.4865 

0 

.2112 

0.6865 

335.1682 

0 

.2162 

0.6797 

334.8472 

0 

.2211 

0.6730 

334.5234 

0 

.2261 

0.6663 

334.1973 

0 

.2310 

0.6598 

333.8685 

0 

.2359 

0.6532 

333.5372 

0 

.2407 

0.6468 

333.2038 

0 

.2455 

0.6404 

332.8679 

0 

.2503 

0.6340 

332.5298 

0 

.2551 

0.6278 

332.1895 

0 

.2598 

0.6216 

331.8472 

0 

.2644 

0.6154 

331. 5027 

0 

.2691 

0.6094 

301.1560 

0 

.  2737 

0.6034 

300,8075 

0 

•  2  7  8  3 

0.5974 

330. 4567 

0 

.  2828 

0.5916 

330.1042 

0 

.  28?3 

0.5858 

329.7497 

0 

,2918 

0.5800 

329.3931 

0 

.2963 

0.5743 

329.0348 

0 

.  3007 

0.5687 

328.6740 

0 

.3051 

0.5632 

328.2875 

0 

.309^ 

0.5577 

327. 9237 

0 

.3137 

0.5552 

327.7587 

0 

.3157 
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WELLBORE 

DEPTH 

(FT) 

0.00 

167.49 

267.04 

366.61 

466.19 

565.78 

665.39 

765.01 

864.64 

964.27 

1063.92 
1163.58 

1263.24 

1362.92 
1462.60 
1562.28 

1661.98 

1761.68 
1861.38 
1961.09 
2060.81 
2160.53 

2260.25 

2359.98 

2459.71 

2559.45 

2659.19 

2758.92 

2858.68 

2958.43 

3058.18 

3157.93 

3257.68 

3357.43 

3457.19 

3556.94 

3656.70 

3756.45 

3856.21 

3955.96 

4055.72 

4155.47 

4255.22 

4354.97 

4454.72 

4554.47 

4654.22 
4753.96 

4853.70 


AFTER  3650.00  DAYS 


steam  temp 

CUM  HEAT 

LOSS 

Quality 

casing  tfmp 

F.  LOcc 

(  F) 

( RTU/HR ) 

(  F) 

467.9000 

O.OOOOE- 

•38 

1.0000 

362.9467 

n.nnnn 

466.0941 

0©  5168E 

05 

0.9865 

362.6183 

0.0089 

465.5442 

0.8222E 

05 

0.9786 

362.4132 

0.0142 

464.9845 

0.1126E 

06 

0.9707 

362.1995 

0.0194 

464.4151 

0.1429E 

06 

0.9629 

361.9776 

0.0246 

463.8362 

0. 1730E 

06 

0.9551 

361.7477 

0.0298 

463.2480 

0®  2030E 

06 

0.9474 

361.5096 

0.0350 

462.6505 

0 ®  2  3 2  8E 

06 

0*9398 

361.2639 

0.0401 

462.0438 

0.2625E 

06 

0.9322 

361.0103 

0.0452 

461,4282 

0.2921F 

06 

0.9246 

360.7492 

0.0503 

460.8036 

0.321 5F 

06 

0.9171 

360.4809 

0.0584 

460.1703 

0®  3508E 

06 

0.9097 

360.205 0 

0.0604 

459.528 2 

0.3799E 

06 

0.9023 

359.9221 

0.0654 

458.8776 

0.4089E 

06 

0.8950 

359.6321 

0.0704 

458.2184 

0.4377E 

06 

0.8877 

359.3350 

0.0754 

457.5509 

0.4664E 

06 

0.8805 

359.0312 

0.0803 

456.8749 

0.4950E 

06 

0.8733 

358.7206 

0.0852 

456.1907 

0.5234E 

06 

0.8662 

358.4034 

0.0901 

455.4983 

0.5516E 

06 

0.8591 

358.0795 

0.0950 

454.7977 

0.5797E 

06 

0.8522 

357.7492 

0.0998 

454.0890 

0.6076E 

06 

0.8452 

357.4124 

0.1046 

453.3722 

0.6354E 

06 

0.8383 

357.0694 

0.1094 

452.6475 

0.6631E 

06 

0.8315 

356.7200 

0.1142 

451.9147 

0.6905E 

06 

0.8248 

356.9644 

0.1189 

451.1740 

0.7179F 

06 

0.8180 

356.0026 

0.1236 

450.4254 

0.7450F 

06 

0.8114 

355.6348 

0.1283 

449.6688 

0.7720E 

06 

0.8048 

355.2609 

0.1329 

448.9042 

0.7989E 

06 

0.7983 

354.8808 

0.1376 

448.1315 

0.8256E 

06 

0.7918 

354.4946 

0.1422 

447.3512 

0.8521E 

06 

0.7854 

354.1026 

0.1467 

446.5629 

0.8785E 

06 

0.7791 

353.7047 

0.1513 

445.7667 

0.9047E 

06 

0.7728 

353.2851 

0.1558 

444.9625 

0.9308E 

06 

0.7665 

352.8751 

0.1603 

444.1502 

0.9567E 

06 

0.7604 

352.4594 

0.1647 

443.3299 

0.9824E 

06 

0.7542 

352.0376 

0.1692 

442.5015 

0. 1008E 

07 

0.7482 

351.6100 

0.1736 

441.6650 

0.1033E 

07 

0.7422 

351.1762 

0.1780 

440.8202 

0..1059E 

07 

0.7363 

350.7365 

0.1823 

439.9671 

0. 1084E 

07 

0.7304 

350.2907 

0.1866 

439.1056 

0.1 109E 

07 

0.7246 

349.8388 

0.1909 

438.2356 

0.1  133E 

07 

0.7188 

349.3807 

0.1952 

437.3570 

0.1158E 

07 

0.7132 

348.9164 

0.1994 

436.4697 

0. 1 182E 

07 

0.7075 

348.4458 

0.2036 

435.5736 

0.1207E 

07 

0.7020 

347.9690 

0.2078 

4  34. 66  8  5 

0.1231E 

07 

0.6965 

347.4855 

0.2119 

433.7542 

0.1255E 

07 

0.6910 

346.9956 

0.2160 

432.8307 

0. 1278E 

07 

0.6857 

346.4990 

0.2201 

431.8977 

0.1302E 

07 

0.6803 

345.9955 

0.2242 

430.9550 

0.1325E 

07 

0.6751 

345.4853 

0.2282 

; 
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LLBORE 

AFTER 

3650.00 

DAYS 

EPTH 

STEAM 

TEMP  CUM 

HEAT 

(  FT) 

(  F) 

( STU/HR 

4953.44 

430.0024 

0.1348E 

5053.17 

429.0398 

0.1371E 

5152.90 

428.0669 

0  ® 1 394  E 

5252.63 

427.0833 

0.1417E 

5352.35 

426.0889 

0. 1440E 

5452.07 

425.0834 

0  ® 1 462E 

5551.78 

424.0664 

0. 1484E 

5651.49 

423.0375 

0.1506E 

5751.19 

421.9966 

0.1528E 

58  50.89 

420,9430 

0, 1549E 

5950.58 

419.8764 

0. 1571 F 

6000.00 

419,3427 

0.1581F 

QUAL I TY 

CASING  TEMP 
(  F) 

F.  LOS 

0.6699 

344,9680 

0.2322 

0.6648 

344,4435 

0.2362 

0.6597 

343.9116 

0.2401 

0.6547 

343.3722 

0.2440 

0.64Q8 

342.8250 

0.2479 

0.6449 

342.2556 

0.2517 

0.6401 

341.6925 

0.2556 

0.6354 

341.1210 

0.2593 

0,6307 

340.5407 

0.2631 

0.6261 

339.9515 

0.2668 

0.6216 

339,3529 

0.2705 

0.6194 

339.0526 

0.2723 

LO 

) 

07 

07 

07 

07 

07 

07 

07 

07 

07 

07 

07 

07 


, 

•  ■ 
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